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ABSTRACT 

We compare line emission calculated from theoretical disk models with optical to submillimeter wavelength 
observational data of the gas disk surrounding TW Hya and infer the spatial distribution of mass in the gas 
disk. The model disk that best matches observations has a gas mass ranging from ~ 10 4 to 10 -5 M 0 for 
0.06 AU < r < 3.5 AU and ~0.06 M 0 for 3.5 AU < r < 200 AU. We find that the inner dust hole (r < 3.5 AU) 
in the disk must be depleted of gas by ~l-2 orders of magnitude compared with the extrapolated surface density 
distribution of the outer disk. Grain growth alone is therefore not a viable explanation for the dust hole. CO 
vibrational emission arises within r ~ 0.5 AU from thermal excitation of gas. [Oi] 6300 A and 5577 A forbidden 
lines and OH mid-infrared emission are mainly due to prompt emission following UV photodissociation of OH 
and water at r <0.1 AU and at r ~ 4 AU. [Neu] emission is consistent with an origin in X-ray heated neutral 
gas at r < 10 AU, and may not require the presence of a significant extreme-ultraviolet (hv > 13.6 eV) flux from 
TW Hya. H 2 pure rotational line emission comes primarily from r ~ I to 30 AU. [0 1 ] 63 /nn, HCO + , and CO pure 
rotational lines all arise from the outer disk at r ~ 30-120 AU. We discuss planet formation and photoevaporation 
as causes for the decrease in surface density of gas and dust inside 4 AU. If a planet is present, our results suggest a 
planet mass ~4-7 Mj situated at ~3 AU. Using our photoevaporation models and the best surface density profile 
match to observations, we estimate a current photoevaporative mass loss rate of 4 x 10 -9 M 0 yr 1 and a remaining 
disk lifetime of ~5 million years. 
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1. INTRODUCTION 

Circumstellar disks are widely believed to go through a tran- 

sition epoch as they evolve from an initial, gas-rich, classical 
T Tauri (CTTS) phase to the weak-lined T Tauri (WTTS) de- 
bris disk phase (e.g., Strom et al. 1989; Skrutskie et al. 1990; 
Najita et al. 2007). Transition disks, which typically are a small 

fraction (~10%-20%) of all disks, are characterized by an 
outer, optically thick dust disk and an inner dust-depleted re- 
gion or “hole.” However, it is unclear if all disks go through 
the transition phase (e.g., Currie et al. 2009; Muzerolle et al. 
2010). Transition disks are sometimes found to be accreting 
at rates similar to young disks, suggesting that the dust holes 
are not devoid of gas, but many objects do not show signs 
of accretion (Sicilia-Aguilar et al. 2010). The duration of the 
transition epoch is also uncertain (e.g., Muzerolle et al. 2010), 
although a recent study of the statistics of spectral energy distri- 
butions (SEDs) suggests times of ~a few 10 5 years (Ercolano & 
Clarke 2010). Studying these partially evacuated disks can po- 
tentially provide many insights into how disks evolve and form 
planets. 

At least three different scenarios have been proposed to 
explain the morphology of transition disks (see Najita et al. 
2007 for a summary). First, the infrared opacity hole may 
not be a true depletion of material, but the result of dust 
evolution. Coagulation and grain growth in disks might lead 
to the formation of larger, “invisible” solid objects, thereby 
emptying the inner regions of small, infrared-emitting dust. In 


this sequence of events, gas and large solids or planetesimals 
remain and there is no true mass depletion (e.g., Dullemond 
& Dominik 2005). In a second situation, the hole in the disk is 
caused by dynamical interactions with an embedded planet (e.g., 
Skrutskie et al. 1990; Marsh & Mahoney 1992; Takeuchi et al. 
1996; Varniere et al. 2006). The Jovian or larger mass planet first 
creates a gap in the disk and then accretes most of the accretion 
flow from the outer disk, while a small fraction continues past the 
planet to accrete onto the star (e.g., Lubow & D’Angelo 2006). 
The third explanation for the inner hole is photoevaporation, 
a process by which high energy radiation from the central star 
heats the gas to escape speeds and creates gaps (and then quickly 
central holes) a few AU from the star (e.g., Alexander et al. 
2006; Gorti et al. 2009; Owen et al. 2010). Determining which 
of the above three mechanisms dominate in the creation of the 
inner opacity hole may be key to deciphering the role played 
by transition disks in disk evolution, planet formation, and disk 
dispersal. 

TW Hya is the closest known (at 51 pc; Mamajek 2005) 
classical T Tauri star with a nearly face-on (e.g., Qi et al. 
2006) transition disk that exhibits an inner hole in its 
(IR-emitting) dust distribution (Calvet et al. 2002; Ratzka et al. 
2007; Hughes et al. 2009). It is a very well-studied object — the 
central star and disk have been investigated from the X-ray 
to centimeter wavelengths. The X-ray, far-ultraviolet (FUV; 
6 eV < hv < 13.6 eV), and optical flux from the central 
star are well known (e.g., Kastner et al. 2002; Herczeg et al. 
2002; Stelzer & Schmitt 2004; Robrade & Schmitt 2006). The 
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dust 6 * * disk has been studied and extensively modeled, with one 
of the earliest detailed models being that of Calvet et al. (2002, 
hereafter C02) who inferred the presence of a dust depleted inner 
cavity of size ~4 AU from the SED. A very small amount of dust 
persists in the inner hole, and the observed 10 jim silicate emis- 
sion feature has been attributed to this component (C02). Ratzka 
et al. (2007) also model the dust disk to infer a dust cavity, but 
with a smaller size, ~0.5-0.8 AU. The presence of a dust hole, 
or at least a significantly dust depleted inner cavity, has been 
confirmed by imaging and interferometric studies (Eisner et al. 
2006; Hughes et al. 2007). The star has a low, but measurable, 
accretion rate with estimates ranging from 4 x 10~ 10 M Q yr~' 
(Muzerolle et al. 2000) to >10~ 9 M Q yr -1 (Alencar & Batalha 
2002; Eisner et al. 2010), indicating the presence of significant 
gas in the inner dust hole. 

The nature of the transition disk around TW Hya is uncertain. 
While a dust hole is clearly present and gas must be present 
in the hole to account for the observed accretion, the amount 
of gas in this region is not well determined. A radial velocity 
planet was reported by Setiawan et al. (2008), but located 
closer to the star (~0.04 AU) than a planet that could cause 
the 4 AU dust edge. (This detection has been disputed and the 
radial velocity variations attributed to stellar surface activity; 
see Huelamo et al. 2008 and Figueira et al. 2010.) TW Hya is 
known to be a strong FUV and X-ray source and observed to 
have outward flows traced by [Ne ii] 12.8 // m emission (Herczeg 
et al. 2007; Pascucci & Sterzik 2009), which implies that 
the disk is presently undergoing photoevaporation. However, 
photoevaporation theory has traditionally predicted a clean 
evacuated inner region devoid of dust and gas soon after 
the formation of a hole, which the observed accretion rate 
onto the central star contraindicates. Owen et al. (2010) argue 
that the inner disk draining epoch, although short (<10 5 years, 
the local viscous timescale at IM- AU where the gap opens), 
may still be observed as an accreting transition disk. 

Clearly, spectral line observations of gas in the TW Hya 
disk might distinguish between the above scenarios. In recent 
years, many detections of gas emission lines from TW Hya 
have become available. The disk around TW Hya was detected 
in the pure rotational transitions of CO by Zuckerman et al. 
(1995), van Zadelhoff et al. (2001), and later by Qi et al. 
(2004, 2006), who also imaged the disk. Herczeg et al. (2002) 
detected UV emission from H 2 from the disk surface near the 
star. Emission in the near-infrared from vibrationally excited H 2 
was detected by Bary et al. (2003). Thi et al. (2004) detected 
submillimeter emission from the molecules CO, HCO + , DCO + , 
CN, and HCN. CO rovibrational emission from the disk has been 
observed (Rettig et al. 2004; Salyk et al. 2007 ; Pontoppidan et al. 
2008). Resolved line emission of [Neu]12.8 pm was reported 
by Herczeg et al. (2007) and Pascucci & Sterzik (2009). Bitner 
et al. (2008), as part of a ground-based survey using TEXES 
(Richter et al. 2003), failed to detect pure rotational H 2 emission 
from the disk, although Najita et al. (2010) detected emission 
with the larger beam of Spitzer. Najita et al. (2010) present a 
high-quality Spitzer Infrared Spectrograph (IRS) spectrum with 
numerous emission features of H 2 , OH, CO 2 , HCO + , [Nen], 
[Ne hi], and H 1 . Thi et al. (2010) recently reported the Herschel 
PACS detection of the [0 1 ] 63 pt m line from TW Hya. In 
addition to this rich data set, strong optical emission has been 
detected from [0 1] 6300 A and 5577 A (Alencar & Batalha 

6 We refer to solid particles with sizes smaller < 1 mm as “dust.” Larger 

solids are ignored as they have too little surface area and are therefore not 

important for gas physics and topics relevant to this paper. 


2002, S. Edwards 2010, private communication). Pascucci et al. 
(2011) reanalyzed the optical data to also report the detection 
of emission from [S 11 ] 4069 A and [Mg 1 ] 4571 A. 

The abundant observational data available for the TW Hya 
system make it an excellent target for studying the physical 
nature of transition disks through theoretical modeling. Gas 
emission line strengths are very sensitive to the density, temper- 
ature, and the abundance of the emitting species (e.g., Gorti & 
Hollenbach 2008). The emission lines observed largely originate 
from trace species, and in order to interpret the data accurately, 
detailed models of gas chemistry and emission are needed. In 
this paper, we use our previous disk models (Gorti & Hollenbach 
2004, 2008; Hollenbach & Gorti 2009; hereafter GH04, GH08 
and HG09, respectively) to compute gas emission line luminosi- 
ties from the disk around TW Hya. In Section 2, we provide a 
short description of the available data and current understanding 
in each case. We briefly describe the modeling procedure and 
then discuss disk models that best match observational data in 
Section 3. We consider the implications of our results on the na- 
ture of the transition disk around TW Hya (Section 4), discuss 
the constraints set on disk mass in Section 5, and finally present 
our conclusions in Section 6. 

2. GAS EMISSION LINES: OBSERVATIONAL DATA 

We summarize below the observed gas line emission at 
various wavelengths. For clarity of discussion, we do not present 
the data in chronological sequence of detection, but instead 
group observations into three wavelength regions, (1) far- 
infrared and submillimeter features, (2) the mid-infrared (MIR) 
features, and (3) the optical and near-infared (NIR) features. 
Table 1 lists the line luminosities for the observed gas emission 
lines. We do not discuss the dust emission observations, but refer 
the reader to C02, Eisner et al. (2006), Ratzka et al. (2007), and 
references therein. 

Sub-millimeter and far-infrared emission. Detection of CO 
2-1 pure rotational emission from TW Hya was first reported in 
the seminal paper on gas disk lifetimes by Zuckerman et al. 
(1995). Several molecular species (CO, HCN, CN, HCO + , 
DCO + ) have been detected in the TW Hya disk at millimeter and 
submillimeter wavelengths using single-dish telescopes (van 
Zadelhoff et al. 2001; van Dishoeck et al. 2003; Thi et al. 
2004). Wilner et al. (2003) presented the first spatially resolved 
observations of the disk in the HCO + J — 1-0 line. This was 
followed by CO rotational line images in the 2-1, 3-2, and 6-5 
transitions using the Submillimeter Array (SMA) by Qi et al. 
(2004, 2006). 

Thi et al. (2010), using the Herschel Space Observatory, 
recently reported the detection of the [0 1 ] 63 p m line and 
set upper limits on the [Oi] 145 /im and [C 11 ] 1 58 /im lines. 
In conjunction with I3 CO measurements, they interpret the 
emission line fluxes as indicative of a low gas/dust mass ratio 
in the disk (~2.6-26) and estimate the disk gas and dust masses 
as (0.5-5) x 1 0 3 M q and 1.9 x 10~ 4 M e , respectively. 

Mid-infrared emission. The [ Ne 11 ] 12.8 pm forbidden line 
was first reported by Ratzka et al. (2007) from Spitzer IRS obser- 
vations along with Hi 6-5 and 7-6 emission lines. The [Nen] 
line was subsequently spectrally resolved using MICHELLE 
on Gemini N by Herczeg et al. (2007) and they found the 
line profile to be consistent with a single Gaussian centered 
at — 2 ± 3 km s' 1 in the stellocentric frame with an intrinsic 
FWHM of ~21 km s -1 . They concluded that the broad line 
width could result from turbulence in a warm disk atmosphere, 
Keplerian rotation of gas located ~0. 1 AU from the star, or 
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Table 1 

Observations of Gas Emission Lines 


Gorti et al. 


Line 

X. 

Line Luminosity (L©) 

Reference 

Obs. Facility 

CO 6-5 

434 /um 

4.6 x 10~ 8 

Qi et al. 2006 

SMA 

CO 3-2 

867 /x m 

2.4 x 10 -8 

Qi et al. 2006 

SMA 

CO 2-1 

1.3 mm 

7.9 x 10" 9 

Qi et al. 2006 

SMA 

HCO+4-3 

841 (tin 

7.7 x 10^ 9 

van Zadelhoff et al. 2001 

JCMT 

Oi 

63 (tm 

2.7 x 10- 6 

Thi et al. 2010 

Herschel PACS 

Oi 

145 (tm 

<4 x 10~ 7 

Thi et al. 2010 

Herschel PACS 

Cii 

157 /xm 

<4.8 x 10“ 7 

Thi et al. 2010 

Herschel PACS 

Neu 

12.8 (tm 

3. 9^1. 7 x 10~ 6 

Najita et al. 2010 

Spitzer IRS 



3.0 x 10- 6 

Pascucci & Sterzik 2009 

VLT/VISIR (r < 10 AU) 

Ne hi 

15.5 /xm 

2.0 x 10" 7 

Najita et al. 2010 

Spitzer IRS 

H 2 S(2) 

12.28 /xm 

4.8 x 10~ 7 

Najita et al. 2010 

Spitzer IRS 



<4.8 x 10~ 7 

Bitner et al. 2008 

Gemini/TEXES (r < 13.7 AU) 

H 2 S(l) 

17.0 /xm 

0.96-1.0 x 10- 6 

Najita et al. 2010 

Spitzer IRS 



<5.5 x 10" 7 

Bitner et al. 2008 

Gemini/TEXES (r < 20 AU) 

OH 

MIR-SH 

6.2 x 10- 6 

Najita et al. 2010 

Spitzer IRS 

COvib 

NIR 

1.0 x 1 0- 5 

Salyk et al. 2007 

Keck II/Echelle 

H2,vib 

2.12 /xm 

7.9 x 10~ 8 

Bary et al. 2003 

KPNO/Phoenix 

H 2 fluo. 

UV 

1.9 x 10- 4 

Herczeg et al. 2002 

HST/ STIS 

Oi 

6300 A 

1.0 x 10~ 5 

S. Edwards 2010, private communication 

Keck/Echelle 

Oi 

5577 A 

1.4 x 10- 6 

S. Edwards 2010, private communication 

Keck/Echelle 

Sii 

4069 A 

3.0 x 10 -6 

Pascucci et al. 201 1 

FEROS 


a photoevaporative flow. These observations were followed by 
high-resolution spectra using VISIR on Melipal/Very Large 
Telescope (VLT) by Pascucci & Sterzik (2009) and more re- 
cently by Pascucci et al. (201 1). They found that the [Ne ii] line 
profile (intrinsic FWHM ~13 km s -1 ) and the ~4- 6 km s' 1 
blueshift with respect to the star are consistent with theoret- 
ical predictions of a photoevaporative flow driven by stellar 
extreme-ultraviolet (EUV) photons (Alexander 2008). The line 
profile may also be consistent with photoevaporation by X-rays 
and the production of [Neu] in predominantly neutral. X-ray 
heated layers (Ercolano & Owen 2010). 

Najita et al. (2010) recently presented high quality Spitzer 
IRS spectra of TW Hya in the 10-20 /rm region. The spectra 
show many emission features including TL S(l) and S(2) pure 
rotational lines, many Hi lines, [Neu], [Nem], and several 
transitions of OH and bands of HCO + and CO 2 . They find no 
strong lines from H 2 O, C 2 H 2 , or HCN as seen in other classical 
T Tauri stars without dust holes (Carr & Najita 2008; Salyk 
et al. 2008; Pascucci et al. 201 1 ; Pontoppidan et al. 2010). They 
conclude that the [Ne 11 ] emission is consistent with origin in 
the X-ray heated disk surface. The detected OH emission arises 
from transitions originating in energy levels as high as 23,000 K 
above ground and they interpret the emission as likely to arise 
from photodissociation of water in the disk since this process 
results in OH in very high rotational states (van Dishoeck & 
Dalgarno 1984a). The relative strengths of the Hi emission are 
found to be consistent with case B recombination in gas possibly 
associated with the accretion flow close to the star. 

Near-infrared, optical, and UV emission. Emission from 
molecular hydrogen in excited vibrational states has been 
detected from TW Hya. Weintraub et al. (2000) detected 
emission from H 2 in the disk around TW Hya in the 1-0 S(l) 
line at 2.1218 p.m through high-resolution spectroscopy using 
CSHELL on NASA’s Infrared Telescope Facility. Herczeg et al. 
(2002) detected H 2 UV emission (with FWHM ~ 1 8 km s -1 ) 
consistent with excitation caused by UV-induced fluorescence 
from a warm surface layer ~2-3 AU from the central star. 
Ingleby et al. (2009) identified a feature in the STIS spectra at 


1600 A which they attribute to electron impact H 2 emission and 
derive a lower limit on the surface density of gas at ~ 1 AU to 
be ~5 x 10~ 5 g cnT 2 . 

CO rovibrational emission from hot gas, presumably in the 
inner (r < 4 AU) disk, has also been detected. Rettig et al. (2004) 
first detected CO from the disk, and concluded that the observed 
emission was optically thin and consistent with a CO gas mass 
of ~6 x 10 21 g. The low emission temperature of ~430 K was 
interpreted as due to a clearing of the inner disk out to a radial 
distance of ~0.5 AU. Salyk et al. (2007) expanded on these 
results with observations that cover a larger range of rotational 
energies. They estimate emission radii for the CO gas to be 
~0.2-l AU, and derive higher temperatures (~800 K) based 
on their excitation diagram analysis. Pontoppidan et al. (2008) 
have recently presented velocity-resolved spectroastrometric 
imaging of the rovibrational lines (4.7 /im) using the CRIRES 
high-resolution infrared spectrometer on the VLT. Keplerian 
disk models are fit to the spectroastrometric position-velocity 
curves to derive geometrical parameters of the molecular gas 
and the CO gas emission from TW Hya is calculated to come 
from r ~ 0.1-1 AU. 

Optical lines from atomic oxygen have been detected from 
TW Hya. Emission from forbidden lines of [Oi] at 6300 A 
and 5577 A have been observed around many young stars with 
disks and outflows (e.g., Hartigan et al. 1995). Apart from a 
high velocity component (HVC, Av ~ 100-200 km s -1 ) that is 
attributed to shocks in the outflows, there is also a low velocity 
component (LVC, Av ~ 10 km s 1 ) that is believed to originate 
in the disk or from the base of a slow disk wind. The origin 
of this line has remained somewhat of a mystery, as physical 
conditions in disks are such that reproducing the observed high 
line luminosities is difficult (e.g., Font et al. 2004; Hollenbach 
& Gorti 2009). The higher excitation [0 1 ] 5577 A line often 
accompanies the LVC, but is never associated with the HVC. The 
optical forbidden lines often exhibit small blueshifts (typically 
~ — 5 ± 2 km s -1 ; Hartigan et al. 1995). Ercolano & Owen 
(2010) explain the LVC emission as arising from an X-ray 
heated, neutral photoevaporative flow that then gives rise to 
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the observed blueshifts. However, the forbidden line emission 
in TW Hya is not blueshifted and does not appear to participate 
in a photoevaporating flow (Pascucci et al. 2011). TW Hya has 
emission lines (only in the LVC) at [Oi] 6300 A and 5577 A, 
with line luminosities estimated at 10~ 5 and 1.4 x 10 6 L 0 , 
respectively (Alencar & Batalha 2002; Pascucci et al. 2011). 
Pascucci et al. (20 1 1 ) obtain an intrinsic FWHM of the [ O i] lines 
as ~ 1 0 km s~ 1 and also report [S ii] 4069 A emission, estimated 
at 3.6 x 10 -6 L q and non-detection of the [S ii] 6731 A line. As 
we show later, very high gas temperatures or densities are needed 
to excite the [Oi] 5577 A line and reproduce the observed ratio 
of the 6300/5577 lines (T gas > 30,000 K for n e < 10 6 cm" 3 or 
Tgas A 4000 K for;;,. > 10 s cm -3 ). Thus, it is a challenge for the 
5577 A emission to arise thermally (however, also see Ercolano 
& Owen 2010). Here, we account for the [O i] emission and line 
ratio as being a result of photodissociation of OH in the disk 
surface layers (e.g., Storzer & Hollenbach 1998; Acke et al. 
2005). 

3. GAS EMISSION LINES: MODELING 

We use our recently developed theoretical disk models to 
model the gas emission from the disk around TW Hya. The 
stellar parameters are reasonably well determined and the dust 
distribution has been previously modeled, and the focus of 
this paper is therefore the gas disk modeling. We vary disk 
parameters to find the model that best fits the observational data, 
with the principal variable being the surface density distribution. 
We find that two disk models with different inner disk masses 
(for r < 4 AU), but with the same outer disk mass (r > 4 AU), 
provide reasonable fits to the data. We describe these disk 
models in what follows and discuss in detail the case with the 
lower inner disk mass, as these results are a somewhat better 
match to the observed line emission. We call this our standard 
disk model (hereafter referred to as SML, “L” for low inner disk 
mass). We also discuss the results from an alternate model with 
a higher mass in the inner disk and refer to that composite model 
as SMH. 

3.1. Model Description and Inputs 

Our theoretical disk models start from a prescribed disk sur- 
face density distribution and stellar spectrum and proceed to 
self-consistently calculate the disk structure from the gas tem- 
perature and density distribution. The gas disk models simulta- 
neously solve for the thermal and chemical properties of the gas 
and vertical hydrostatic equilibrium. Gas and dust temperatures 
are computed separately and can deviate significantly in the sur- 
face layers where most of the line emission originates. The gas 
radiative transfer follows an escape probability formalism. Dust 
radiative transfer is done using a modified two-layer approach 
(e.g., Chiang & Goldreich 1997; Dullemond et al. 2001 ; Rafikov 
& De Colie 2006). Our chemistry consists of ~85 species and 
~600 reactions, includes photodissociation and photoionization 
by EUV, FUV, and X-rays, and is focused toward accuracy in the 
disk surface layers. Our standard models do not include freeze- 
out of species on grain surfaces, but we test for the effects of 
freeze-out where relevant below. We also assume that H 2 forms 
on grain surfaces in our disk models for all grain temperatures 
below ~400 K (Cazaux & Tielens 2004). For further details, 
see GH04, GH08, and HG09. 

Stellar parameters and the incident radiation field are key 
inputs to our theoretical models and are reasonably well de- 
termined for TW Hya. The source is at a distance of 51 pc 



Figure 1 . Input spectrum used for TW Hya, which is a composite of the stellar, 
FUV, and X-ray spectra shown here as the flux at the stellar surface as a function 
of photon energy. The vertical dashed lines indicate the EUV region which has 
been ignored in our modeling. 

(Mamajek 2005) and the disk is almost face-on with a very low 
inclination angle (7° for the outer disk, Qi et al. 2004; 4° for the 
inner disk, Pontoppidan et al. 2008). We adopt M = 0.7 M 0 , 
R = 1 R q , and T = 4200 K for the stellar mass, radius, 
and effective temperature (Webb et al. 1999). We construct 
a composite spectrum for the star by adding observed X-ray 
and FUV spectra. We use the Far-Ultraviolet Spectroscopic 
Explorer / International Ultraviolet Explorer spectrum from the 
NASA HEASARC archive (Valenti et al. 2003; Herczeg et al. 
2002). The FUV region is split into several energy bins for ac- 
curately computing gas opacity as described in GH04. Herczeg 
et al. (2002) derive a reconstructed Lya flux ~30% higher than 
the observed flux and we use this value for the Lya band of 
our FUV spectrum. The total FUV (6-13.6 eV) flux is ap- 
proximately 3 x 10 31 erg s _1 . The X-ray spectrum is from 
the XMM-Newton archive (Robrade & Schmitt 2006) and ex- 
tends from 0.3 to 5keV, and with a total X-ray luminosity of 
1.3 x 10 30 erg s 1 . We extrapolate the observed spectrum to 
lower and higher energies (0.1-10 keV) and this results in a 
total X-ray luminosity of 1.6 x 10 30 erg s'. Although the star 
might also have an EUV component (13.6-100 eV), perhaps in- 
dicated by the observed photoevaporative flow (Alexander 2008; 
Pascucci & Sterzik 2009), its strength can only be crudely es- 
timated by indirect arguments (Alexander et al. 2005) and we 
do not include this uncertain contribution to the high energy 
flux. We discuss below, where relevant, the effect on line emis- 
sion fluxes if an EUV flux were to be added. Figure 1 shows the 
composite spectrum used as model input and the vertical dashed 
lines in the figure indicate the neglect of an EUV radiation field. 

Two different dust disk models have been proposed for 
TW Hya. by Calvet et al. (2002) and Ratzka et al. (2007). 
C02 developed detailed models of the TW Hya dust SED to 
constrain many aspects of the disk structure, and these models 
have been validated by subsequent observational studies (e.g., 
Qi et al. 2004; Eisner et al. 2006; Hughes et al. 2007). The flux 
deficit at <20 // m and a sharp increase or excess at "-20-60 /im 
led them to propose the presence of an inner hole in the disk 
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Table 2 

Dust Model Parameters 


Optically Thin Inner Disk 


Inner radius 

0.06 AU 

Outer radius 

3.5 AU 

Dust disk mass 

2.4 x 10 25 g 

Surface Density X 

a r~ l 

Minimum grain size 

0.9 /zm 

Maximum grain size 

2.0 fim 

Dust composition 

Glassy Pyroxene 

Optically Thick Outer Disk 

Inner radius 

3.5 AU 

Outer radius 

200 AU 

Dust disk mass (with a < 1 mm) 

4.8 x 10 29 g 

Surface Density X 

Composite (Figure 1) 

Minimum grain size («min) 

0.01 fun 

Maximum grain size (a max ) 

1 mm 

Dust composition 

Silicates + graphite 


at r < 4 AU. The deficit is ascribed to a lack of dust in the 
hole, and the excess arises from direct illumination of the rim 
of the optically thick dust disk beyond 4 AU. There is also a 
small mass of dust, 0.5 lunar masses, in the inner cavity which 
gives rise to the observed spectral feature at 10 /im. Ratzka et al. 
(2007) situate the inner edge of the outer optically thick dust 
closer to the star at ~0.5-0.8 AU. They model both the SED and 
the 10 fim visibilities in detail by including sedimentation in the 
disk, where the dust is reduced in the upper atmosphere due to 
settling processes. Our dust models at present do not include 
settling. We defer gas modeling of the disk using the Ratzka 
et al. (2007) dust distribution to future work, and keep our dust 
model relatively uncomplicated at present. 

We avail ourselves of the dust models of C02. For simplicity, 
we adopt many features of this dust model, specifically the 
dust size and spatial distribution in our standard dust model. 
We consider their favored model with a maximum grain size 
a max = 1 cm and disk mass 0.06 Mq, which corresponds to 
a total dust mass of 7.8 x 10 ~ 4 M Q for their assumed gas/ 
dust ratio. As we ignore dust particles with sizes larger than 
1 mm which are not very relevant for gas physics, our model 
disk “dust” ( m ax = 1 mm) mass is reduced accordingly to 
2.4 x 1CT 4 M 0 (for a grain size distribution n(a ) oc a~ 3 5 , 
where a is the grain size). Our chemical composition of dust 
for the outer disk is simpler than the detailed models of C02. 
We use a mixture of silicates and graphite for the optically 
thick disk and verify the assumption by the SED fit from the 
resulting model dust continuum emission. We note that the 
disk structure (especially in the regions where gas emission 
lines originate) is eventually determined by the gas temperature 
structure, although the dust distribution is sometimes important 
for the heating and cooling of gas. We use pyroxene grains 
for the inner dust hole as suggested by C02 to better fit the 
continuum SED and the 10 /im feature. As we show later, the 
amount of dust present in this inner region is in any case too 
low to significantly affect the gas thermal balance or chemical 
calculations. The details of the adopted dust model are listed in 
Table 2. 

Polycyclic aromatic hydrocarbons (PAHs) can efficiently heat 
gas via the grain photoelectric heating mechanism (Bakes & 
Tielens 1994) and gas heating depends on the assumed PAH 
abundance. Disks around low mass stars are believed to be 
deficient in PAHs by a factor of ~ 10-100 (Geers et al. 2006). 
We typically scale the PAH abundance with the decrease in 


cross-sectional area of the dust per H atom compared with the 
interstellar value. We adopt a PAH abundance of ~ 10~ 7 per H in 
the interstellar medium (e.g., Habart et al. 2003). Therefore, the 
abundance of PAHs in the optically thick disk (r > 4 AU) around 
TW Hya is ~ 10~ 9 . Such a depleted abundance means that X-ray 
heating usually dominates FUV-induced PAH heating of gas in 
the outer disk. The abundance of small grains (a <0.1 /im) in 
the inner disk (r < 4 AU) is either negligible (C02) or very small 
(Eisner et al. 2006). It is therefore likely that the abundance of 
very small dust like the PAHs is low in the inner disk. Such low 
PAH abundances do not influence the heating in the inner disk 
and the PAH abundance in the inner disk is therefore set to zero. 

We do not consider turbulence in the disk beyond 4 AU and all 
cooling line widths are thermal in the optically thick outer disk. 
This assumption is supported at large radii by submillimeter 
observations and modeling of molecular emission from TW Hya 
by Qi et al. (2006) and Hughes et al. (2011). Hughes et al. (201 1) 
infer a small turbulent linewidth that is approximately a tenth 
of the thermal linewidth in the outer disk. The situation may 
be different in the inner disk, r < 1-2 AU. Salyk et al. (2007) 
assume higher turbulent linewidths (v tur b ~ 0.05t >k ~ 4kms _1 , 
where the Keplerian velocity Vk ~ 80 km s _ 1 at 0. 1 AU) in their 
analysis of the rovibrational CO lines from the disk. CO rov ib 
emission arises from hot gas, and the thermal speed of CO for 
~2000 K gas is ~1 km s . The expected thermal linewidth 
is therefore ~2 km s'. Observed line widths are typically 
~15 km s' 1 (Salyk et al. 2007), with presumably a large 
Keplerian contribution for gas located at small radii (~0. 1 AU) 
in a face-on disk with a small inclination (4° ; Pontoppidan et al. 
2008). We follow Salyk et al. (2007) in assuming i> tur b = 0.05i>^ 
in the inner disk (r < 4 AU) regions. 

We also add viscous heating as a mechanism to heat gas in 
the inner disk. The dust-depleted inner cavity might still contain 
dense gas and viscous heating may be an important contributor 
to heating the dust-poor gas, especially in the shielded interior 
near the midplane. We follow the prescription given in Glassgold 
et al. (2004, Equation (12) of that paper) for the viscous heating 
term. The viscosity parameter a is unknown, but is set in our 
models by assuming steady state accretion onto the star at the 
observed rate of ~ 1 CT 9 M e yr" 1 and from the solutions obtained 
for the density and temperature of the gas. 

The surface density distribution of gas in the disk is the 
primary input to our disk models. We begin with the dust 
distribution of C02 which has a rapid increase in surface 
density over ~3.5-4 AU. We then vary the gas surface density 
distribution, Xtr), so that the resulting gas emission lines 
provide a good match to observational data. The resulting 
profile for the SML disk is shown in Figure 2. The SML disk 
has a gas mass ~10 ~ 5 M 0 within 3.5 AU and 0.06 M e for 
3.5 AU < r < 200 AU. Although the dust surface density drops 
by ~10 3 inward of r ~ 4 AU, the gas surface density drops 
only by ~ 1 00 for the SML disk. The alternate model SMH 
has a gas mass ~10 4 M 0 within 3.5 AU, and in this case the 
gas surface density only drops by ~10. The outer disk surface 
density distribution is the same for both SML and SMH cases. 
We discuss later in Section 3.3 how deviations from the models 
SML and SMH worsen the agreement between model results 
and observational data. 

For descriptive purposes, we divide the disk into three distinct 
regions, (1) the inner disk (r < 3.5 AU), (2) the mid-disk region 
(3.5 AU < r < 20 AU), and the (3) outer disk (r > 20 AU). 
Our modeling procedure first constrains the gas in the inner disk 
region, as this gas provides opacity to stellar radiation, shields 
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Figure 2. Surface density distribution of the SML (solid line) and SMH (dashed 
line) disks. The surface density follows a power-law (Jar 1 ) within 4 AU, rises 
over a “rim” region (3.5—4 AU), and is approximately a power-law (X a r °' 7 ) 
from 4 40 AU. In the outer disk (r > 40 AU), the surface density rapidly 
decreases with radius. See Section 3 . 1 for the analytic formulation of X(r). 


gas in the partially exposed rim region at 3.5-4 AU, and affects 
the structure of the outer disk. The inner region of the SML 
disk is not only depleted in dust but also in gas, by >2 orders 
of magnitude compared to that expected if the outer optically 
thick disk was smoothly extrapolated radially inward. A power 
law is assumed for the inner disk surface density distribution, 
l,(r) oc r” 1 ; 0.06 AU < 3.5 AU. We place the inner disk 
edge at 0.06 AU, at the location of the magnetospheric truncation 
radius (Eisner et al. 2006). As there is continued accretion onto 
the star, presumably by funnels along magnetic field lines, this is 
a logical choice for this parameter. Between 3.5 and 4 AU, XO) 
rises sharply and these radii constitute the “rim” region. The 


functional form for £(r) is modeled to fit the appearance of the 
rim of a photoevaporating disk (Gorti et al. 2009) and is given 
by logZ(r) = 18 — exp(ll/r A u) g cm -2 ; 3.5 < r AU < 4. Here, 
r AU = r/1 AU. The dust disk turns optically thick over this 
region. Beyond 4 AU, the surface density is assumed to follow a 
self-similar density profile, as expected from a photoevaporating 
and viscously evolving disk (e.g., Hughes et al. 2007; Gorti et al. 
2009). For r > 4 AU, Z(r) = 500r A °' 7 exp(-r^/100) g cm" 2 . 

3.2. Model Results: Disk Structure 

We briefly discuss the disk density and temperature structure, 
the dominant heating and cooling mechanisms at three represen- 
tative disk radial positions, and then describe the line emission 
calculated from the disk models. 

The density and temperature distribution as a function of 
spatial location in the SML disk is shown in Figure 3. Gas tem- 
peratures in the surface regions where emission lines originate 
are typically ~2000-200 K in the inner and mid-disk regions, 
r < 30 AU. The near absence of dust in the inner disk and the 
sudden rise in dust surface density at the rim causes an increase 
in gas heating and results in a small vertically extended region 
at this radius ~4 AU. However, gas in the inner disk provides 
opacity that shields the midplane regions at ~4 AU, and there 
is only a slight “puffing-up” of the inner rim. (We find that for 
the upper layers of the outer disk where the gas temperature 
decouples from the dust temperature, the column density of gas 
to the star at a given angle (height of the rim) changes by less 
than 10% as the mass of the inner disk changes from 10~ 5 M Q 
to 10 -4 Mq.) As seen in Figure 3 (see Ay = 1 surface), the disk 
structure is overall quite flat, with only modest flaring of the 
outer disk. Gas at depths where Ay ~ 0.1-1 to the star, and at 
radii r > 50 AU, only attains temperatures of ~50-200 K and 
the density in these warm layers is low, n ~ 10 5 — 10 7 cm -3 . 

Disk gas temperature and density is shown as a function of 
height at typical positions in the three different disk regions, 
inner disk at r = 0.1 AU, mid-disk at r = 6 AU, and the outer 
disk at r = 30 AU in Figures 4-6. 

In the inner disk (Figure 4), heating is mainly by X-rays 
all through the disk, and in the midplane FUV is significant 


log n log T 



0 20 40 60 80 100 0 20 40 60 80 100 

Radius (AU) Radius (AU) 


Figure 3. Density and temperature structure of the disk for the SML disk. The outer disk (r > 40 AU) density structure shows a somewhat flattened disk structure, 
resulting in relatively low disk temperatures in these regions. The black line in both panels shows the height at which the visual extinction to the star Ay = 1 . The 
insets show the structure in the inner disk for r < 1 AU where most of the line emission is produced. 

(A color version of this figure is available in the online journal.) 
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Figure 4. Vertical temperature (solid line) and density (dashed line) distribution 
is shown at a typical radius in the inner disk (0.1 AU). The dominant heating 
(above) and cooling (below) mechanisms are shown and the radial extent where 
they are important in overall thermal balance is indicated. The photodissociation 
fronts of CO and H 2 are also marked on the plots. 


only via chemical heating (primarily due to photodissociation 
of water and ionization of carbon). However, to attain the high 
surface temperatures achieved in this region, FUV photons are 
essential to photodissociate molecules capable of efficiently 
cooling the surface gas. The gas temperature increases when H 2 
photodissociates and molecular (H 2 O and OH, formed via H 2 ) 
cooling decreases, and also when CO photodissociates and can 
no longer cool the gas efficiently. In the inner disk, the transition 
from C to CO occurs higher up than the transition from H to H 2 , 
as is typical in low dust regions (e.g., Glassgold et al. 1997). The 
inversion in height can be seen in Figure 4 where the H 2 and 
CO photodissociation fronts are marked. Rovibrational emission 
of CO (CO r ovib) dominates the cooling with some additional 
cooling by H 2 O and OH, until CO photodissociates. At higher 
z, 0 1 6300 A and Lya are the main coolants. 

Figure 5 shows the disk vertical structure in the mid-disk 
at r — 6 AU. In these regions where dust is abundant, dust 
collisions keep gas at the ambient dust temperature in the 
optically thick midplane region. Near the optically thin surface 
layer and near the H 2 photodissociation front, heat released 
from the formation of H 2 is significant. Chemical heating due 
to photoreactions (FUV photodissociation of water) and FUV- 
induced grain photoelectric heating by PAHs is of significance 
over a limited column (Ay ~ 0.1-1). At the surface. X-rays 
dominate the heating. Cooling is mainly due to CO and H 2 in 
the molecular layer, and Oi 6300 A, [Oi] 63 /mi, [Neil], and 
[Arn] fine structure lines cool the surface gas. 

The outer disk vertical structure at r = 30 AU is shown in 
Figure 6. Dust is thermally coupled to the gas by collisions in the 
midplane where gas and dust temperatures are equal. Near the 
Ay — 1 layer, heating by PAHs, photodissociation reactions and 
formation and dissociation of H 2 are important. X-rays dominate 
the heating at vertical column densities of < 1 0 22 cnT 2 . FUV 
photoionization of carbon contributes to the heating (~15%) in 
the region above the CO photodissociation layer. Main coolants 
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Figure 5. Disk structure in the mid-disk region at a radius of 6 AU. Here, the 
dominant heating mechanism is FUV heating by PAHs and photoreactions near 
the A v = 1 layer, and X-rays at the surface. Near the midplane, dust collisions 
are important. The main gas coolants are CO and H 2 rotational lines deeper in 
the disk, and 0 1 (63 /cm and 6300 A ) and [Ne n] and [Ar 11 ] at the disk surface. 
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Figure 6. Disk structure at a representative radius in the outer disk (30 AU). The 
dominant heating and cooling mechanisms and their radial extent is indicated 
in the figure. 

in the outer disk are CO in the molecular region and [0 1 ] 63 /tm 
fine structure emission at the surface. 

3.3. Model Results: Line Emission 

We adopt the following procedure to determine the surface 
density distribution that is described in Section 3.1. In order 
to obtain a gas disk surface density model that best matches 
observed gas line emission, we first constrain the surface density 
distribution in the inner disk, i.e., within 4 AU. This is because 
gas in the inner disk can shield the outer disk regions, affecting 
the density and temperature structure of the disk beyond 4 AU, 
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Table 3 

Model Calculations for Lines (Sum of Inner and Outer Disk Contributions) that Constrain the Inner Disk Mass 


Line 

Obs.Data 

SMH 

SML 

§ 

11 

0 

1 

1 

Contribution 



§ 

11 

0 

1 

1 

(Mi = 10“ 5 M a ) 


from r > 4 AU 


Uo) 

<i 0 ) 

(£0) 

a*©) 

(£©) 

COvib (P Series) 

10~ 5 

2.1 x 10" 5 

1.8 x 10“ 5 

1.6 x 10“ 5 


OH (prompt) 

1.5 x 10" 5 

7.8 x 10“ 6 

4 x 10" 6 

2.8 x 10” 6 

2.8 x 10 -6 

0 1 6300 A 

10~ 5 

1.5 x 10" 5 

8.4 x 10” 6 

8.9 x 10" 6 

6.1 x 10 -6 

( Thermal ) 


(3.1 x 10 -6 ) 

(2.3 x 10“ 6 ) 

(2.8 x 10” 6 ) 

(1.1 x 10“ 6 ) 

(Prompt) 


(1.2 x 10~ 5 ) 

(6.1 x 10~ 6 ) 

(6.1 x 10 -6 ) 

(5 x 10“ 6 ) 

0 1 5577 A 

1.4 x 10~ 6 

1.92 x 10~ 6 

1.6 x 10~ 6 

9.4 x 10~ 7 

7 x 10- 7 

(Thermal) 


(2.2 x 10“ 7 ) 

(1.0 x 10“ 7 ) 

(1.3 x 10 -7 ) 

(1.0 x lO^ 10 ) 

( Prompt ) 


(1.7 x 10~ 6 ) 

(1.5 x 10“ 6 ) 

(8.1 x 10" 7 ) 

(7 x 10~ 7 ) 

H 2 S(2) 

4.8 x 10- 7 

6.8 x 10” 7 

3.4 x 10” 7 

2.5 x 10“ 7 

2.2 x 10 -7 

h 2 S(l) 

1.0 x 10“ 6 

1.3 x 10“ 6 

1.2 x 10~ 6 

1.1 x 10~ 6 

1.1 x 10" 6 

Ne 11 

5.0 x 10“ 6 

3.1 x 10~ 6 

3.0 x 10~ 6 

3.4 x 10~ 6 

2.3 x 10" 6 

Ne hi 

2.0 x 10 -7 

3.0 x 10" 7 

3.0 x 10“ 7 

3.5 x 10 -7 

1.0 x 10“ 7 

O i linewidth (km s -1 ) 

10.0 

11.0 

6.0 

6.1 


a (vise, parameter) 


0.01-0.02 

0. 1-0.2 

1-2 



where the dust disk is optically thick. There clearly is some 
gas in the inner hole, based on the observed accretion onto the 
central star, and on observations of CO rovibrational lines and 
Hj ultraviolet lines. We discuss these emission line constraints 
and the undetected water lines below, and their implications for 
the gas mass in the inner hole. 

Gas opacity. Although the dust is optically thin to stellar ra- 
diation, the column density of gas may be sufficient to provide 
significant optical depth in the inner disk. In fact, the gas may 
be optically thick to its own infrared emission. We use our disk 
chemical models to obtain chemical abundances for gas in the 
inner disk and follow Tsuji (1966) in calculating the monochro- 
matic opacity due to each gas species. At the likely pressures 
and temperatures of our inner disks, H - and H 2 bound-free and 
free-free, CO and H 2 O all can provide significant opacity. H~ 
and Hj" bound-free and free-free continuum opacity can be sig- 
nificant at wavelengths shorter than about 4 pm, and produce 
near IR radiation from the inner edge of the disk at ~0.06 AU. 
The ground state CO vibrational band occurs near 4.7 pm and 
the blend of CO rovibrational lines provide opacity near 5 pm. 
Beyond about 6 pm, the blend of pure rotational H 2 O lines 
dominate the opacity. 

Later on in this subsection, we discuss CO and LLO lines in 
some detail. Here, we make more general arguments based on 
the expected continuum emission from gas in the disk. We use 
the opacities from Tsuji (1966) to show that if the surface density 
of the outer disk is extrapolated inward to produce a ~ 10 -3 M G 
inner disk (we call this the “full disk” or the “undepleted disk”), 
observable emission would result in the wavelength region 
5-15 pm due to the CO and H 2 O emission. We find, by a similar 
analysis, that the H and Hj" bound-free and free-free produce 
excess emission at shorter wavelengths (~2-3 / um) where the 
stellar emission is significant. We use the specific opacity (per 
molecule) for each opacity contributor from Tsuji (1966) and 
calculate the opacity at each radial annulus by multiplying this 
with the vertical column density as obtained from our model disk 
chemical calculations. For illustration, we calculate the flux at 
one wavelength and compare that with the observed emission. 
We choose a wavelength of 8 pm where the stellar contribution 
is negligible, there is no contamination from the silicate feature 
and where there is a significant dip in the continuum emission 
(see C02). At each radial annulus, the photosphere (at 8 /im) is 
defined to lie at the height (z) above the midplane where the 


total optical depth at 8 pm becomes of order unity. Near the 
inner edge of the disk for 0.06 AU < r < 0.08 AU, H and 
Hj" free-free are the main opacity contributors. This opacity 
drops rapidly with decreasing electron density further out in 
the disk. At 0.08 AU < r < 0.2 AU, H 2 O lines dominate the 
8 pm opacity. Having defined the photosphere of the disk from 
the gas opacities, we estimate the 8 pm flux from this surface, 
where the gas temperature ranges from ~300 to 2000 K. The 
disk becomes optically thin to the midplane (X ~ 8 pm) at 
~0.2-0.3 AU. We obtain an 8 pm flux of ~1 Jy, a factor of ~4 
higher than observed. A full gas disk is therefore incompatible 
with the SED. The H 2 O lines at 8 pm therefore yield an estimate 
of the upper limit on gas mass as 2.5 x 1 0 4 M 0 . We discuss 
later our model results for H 2 O line emission and compare them 
to Spitzer IRS upper limits to the line fluxes to obtain a better 
limit on the inner disk mass. 

A similar opacity analysis at a wavelength of 5 pm suggests 
that CO rovibrational emission in excess of what is observed will 
be produced by a full gas disk. We use our model calculations 
and compare them with observed rovibrational CO line strengths 
near 5 pm to estimate the actual mass of the inner disk to be 
~10 -4 to 1CT 5 M 0 . 

In what follows, rather than using an opacity argument, we 
use our disk model results on emergent line fluxes from the 
disk surface to infer the disk mass. We describe the emitting 
regions for each gas emission line diagnostic and also discuss 
the constraints set in each case, where applicable. The outer 
disk gas is attenuated by the gas in the inner disk as described 
earlier. We calculate the outer disk structure using the surface 
density distribution of the inner disk in the SML model and 
use the same results for the outer disk component of the SMH 
model as well, as the attenuation columns of the gas provided at 
4 AU are similar in both of these cases (within ~10%). As our 
disk models are computationally intensive, we do not calculate 
disk models for the outer disk in the SMH case separately. 
In the rest of the paper, the emission in the mid and outer 
disk regions refer to that calculated using gas optical depths 
and absorption columns for the 10~ 5 M Q inner disk mass case 
(SML). Tables 3 and 4 summarize the contribution of the inner 
and outer disk to the emission lines for the SML and SMH 
model disks. Figure 7 shows the radial extent of each of the gas 
emission lines discussed below. The cumulative line luminosity, 
L(r), integrated up to the radius r is shown for each emission 
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Figure 7. Cumulative line luminosities (L(r)) normalized to their total values 
for the SML as a function of disk radius. The molecular lines, including prompt 
emission, usually arise from near Ay = 1 layer, whereas the Oi and Neil 
thermal emission originate above this layer in the disk surface. The 10% and 
90% levels (dotted lines) mark the radii between which most of the emission 
originates. For example, the [Ne n] 12.8 /im line mainly arises from gas located 
between 1 and 8 AU. 

(A color version of this figure is available in the online journal.) 


Table 4 

Calculated Model Line Luminosities for Lines Originating Beyond 4 AU 


Line (X) 

Obs. Data 
(Lq) 

Model 

Uo) 

Oi(63 (im) 

2.7 x nr 6 

3-6 x 10- 6 

Oi(145 (jm) 

<4 x ur 7 

1-3 x 10" 7 

C ii(158 (im) 

< 4.5 x ur 6 

3 x 10" 6 

CO 6-5(430 /im) 

4.6 x ur 8 

8 x 10" 8 

CO 3-2(866 /im) 

2.4 x ur 8 

2.4 x 10 -8 

CO 2-1 (1.3 mm) 

7.9 x ur 9 

7.1 x 1(T 9 

HCO + 4~3(841 Atm) 

7.7 x ur 9 

4 x 1(T 9 


line and is normalized to the total calculated value (from inner, 
mid, and outer disk regions). Also marked are the 10% and 
90% luminosity levels. We quote the corresponding radii below 
for each emission line and this defines the radial extent of the 
emission. 

CO rovibrational emission. CO rov ib emission typically arises 
from a height in the inner disk where the CO is just beginning 
to photodissociate. Gas temperatures in the disk generally 
increase with height and gas temperatures in the emitting 
layers are high ( 7’„ as > 500 K). In dust-poor gas, the CO and 
H 2 photodissociation layers are inverted, as first reported by 
Glassgold et al. (2004), that is, the carbon becomes CO at larger 
z than the H becomes H 2 . We similarly find that the H 2 is 
photodissociated to deeper layers and that CO is coexistent 
with large columns of atomic hydrogen, creating conditions 
favorable for exciting CO rov jb emission (Glassgold et al. 2004). 
The primary formation route to CO is the oxygenation of 
atomic carbon (C+OH). In the layers shielded from FUV 
photodissociation CO is destroyed by collisions with helium 
ions. 


CO vibrational lines are an important constraint for setting 
the gas mass in the inner disk. A high disk mass results in 
high total CO vibrational emission if the lines are optically 
thick, as the luminosity in these lines increases with temperature 
which increases (albeit weakly) with increasing gas mass. The 
total luminosity of the observed lines (select P and R branch 
transitions) is constrained by observations to be <10~ 5 L Q 
(Salyk et al. 2007). We find that inner disk masses of 10~ 4 , 1 O ' 5 , 
and 1 0 6 Mq all give total CO luminosities within a factor of 
~2 of the observations and are all hence viable choices. The CO 
emission is found to be mainly a surface effect dominated by 
hot gas at the surface of the disk, resulting in approximately 
the same total CO luminosity for all three disk masses. To 
distinguish between these cases, the line luminosity of each 
of the CO rovibrational lines needs to be compared with data. 
We estimate the CO rovibrational line luminosities of the v = 
1-0 P(J) series as follows. We adopt the analysis of Scoville et al. 
(1980) and perform an approximate non-LTE calculation for the 
populations of the CO rovibrational levels. We only consider the 
v = 0 and v = 1 states as higher vibrational states are unlikely 
to be significantly occupied at typical disk temperatures. We 
assume that collisions result in a change in the vibrational 
state, and vibrational level populations are determined from 
detailed balance assuming an “effective” Einstein A-value from 
v = 0-1 as described in Najita et al. (1996). We then assume 
that the rotational levels within each v state are thermally 
populated (Scoville et al. 1980). Collisional rates with e~ , H, 
and H 2 are taken from Najita et al. (1996). Oscillator strengths 
and transition frequencies are from the compilations of Kirby- 
Docken & Liu (1978). The optical depth for each transition is 
calculated as detailed in Hollenbach & McKee (1979) from the 
column density of CO molecules in a given rovibrational state, 
and an escape probability formalism is used for line transfer. 

Salyk et al. (2007) detect 15 transitions covering rovibrational 
energies E u /k from 3000 K to 6000 K, where E u is the energy 
of the excited upper level. We also compare our results with 
earlier CO rovibrational detections by Rettig et al. (2004) in 
complementary transitions. There was a typographical error 
in the Rettig et al. (2004) table of line fluxes, and the actual 
fluxes are a factor of 10 below what has been reported. We 
also do not include the P(21-25) transitions of that paper, 
as these data are prone to errors due to photospheric and 
telluric corrections (S. Brittain 2011, private communication). 
We compare our model results with the data in Figure 8. 
We find that a gas depletion from the full disk of 0.1-0.01 
corresponding to a gas mass of 1 0 4 to 10~ 5 M G matches the 
observed line emission fairly well, with the lower mass inner 
disk model (SML) being a better fit to the CO data. A factor 
of 10 lower in mass results in CO emission from warmer gas 
that produces excess emission at high J levels. All of the CO gas 
emission comes from r < 0.5 AU (due to rapidly decreasing gas 
temperatures beyond), in agreement with both the analysis of 
Salyk et al. (2007) and the spectroastrometric imaging study by 
Pontoppidan et al. (2008). Our gas mass estimate for the SML 
disk is consistent with the surface density estimate of Salyk 
et al. (2007) of 1 g cm -2 at r ~ 1 AU where they used a slab 
model, although in our models the emission is mainly from hot 
gas in the disk surface layers. Salyk et al. (2007) also assume an 
H 2 /CO conversion ratio of 5000 as observed toward a massive 
protostar, and caution that their derived mass estimate is likely be 
a lower limit. In our disk models, the hydrogen is predominantly 
atomic in the regions where CO rov ib emission originates (also 
see Glassgold et al. 2004). We explicitly calculate disk chemical 
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Figure 8. Calculated CO rovibrational fluxes for the v = 1-0, P (J = 1-30) 
transitions and the SML inner disk mass ( solid line) are shown for each rotational 
quantum number J. The filled circles are the photosphere and veiling corrected 
measurements from Salyk et al. (2007). Open squares are corrected data from 
Rettig et al. (2004; see the text). The dashed line is for a higher disk mass 
(10 -4 Mq) that also fits the data reasonably well. The 1 0 6 Mr. inner disk mass 
model (dotted line) is a poor fit to the data. 

abundances and find that the H/CO ratio in the hot inner disk 
regions where CO v ib emission arises is ~7000. The gas/dust 
mass ratio in the inner disk is ~500-5000 for the dust mass of 
the Calvet et al. (2002) disk, indicating a relative depletion of 
the dust with respect to the gas at r < 4 AU. 

HiO rotational emission. Our gas opacity analysis above 
suggests that a full undepleted disk would produce excess 
H2O rotational emission in lines whose wavelengths range 
from 6 to 15 pm compared with observations. As we specif- 
ically calculate water abundances and the disk density and 
temperature structure, we can compare the model disk ThO 
rotational line fluxes to the corresponding upper limits from 
Spitzer IRS in the 1 0- 1 9 p m wavelength region (Najita et al. 
2010). No water emission has been detected and line lumi- 
nosities are therefore expected to be lower than ~10 - 7 l 0 . 
We find that for the full undepleted inner disk with a mass 
of 1 0 3 Mq, many rotational lines in the 10-19 p m region are 
above the detection limit. Some of the strongest transitions are 
7 6 i-6i 6 , 762-6i5, 854— 7 o7 , 853-726, and 844-717, all with lumi- 
nosities in the range of (1-3) x 1 0 6 L 0 , compared with the 
upper limits set by Spitzer observations of ~ 10 - 7 ^o . From the 
water rotational line emission constraint alone, we can set a limit 
on the inner disk mass as ~ 1 0 4 Mq, at which the water emis- 
sion is barely detectable. This limit is consistent with the upper 
limit set by CO emission in the preceding discussion, where we 
fix the inner disk mass to be in the range 1 0 4 to 10~ 5 Mq based 
on CO rovibrational emission. 

Hi fluorescent emission. We do not model FF UV fluores- 
cence but draw on the results of Herczeg et al. (2004), who 
conducted a very detailed analysis and modeled the observed 
emission. Our SML disk corresponds to their thick disk model, 
and most of the emission in this case arises from the illumi- 
nated edge of the disk at r — 0.06 AU. Herczeg et al. (2004) 
estimate a hot ( 7„ as ~ 2500 K), H2 mass of ~10 19 g from 


their analysis. We calculate the mass of H2 gas at these tem- 
peratures in our inner disk to be ~3 x 10 19 g, in reasonable 
agreement with the Herczeg et al. (2004) results. All of this gas 
lies very close to the inner edge at ~0.06 AU. At this radius the 
Keplerian speed of the gas from the almost face-on disk ( i =4°) 
is about ~7 km s -1 , the thermal speed is ~4 km s -1 , bringing 
the expected FWHM to ~16 km s _1 very close to the observed 
FWHM of 18 km s -1 (Herczeg et al. 2002). We note that the 
inner disk structure near the magnetospheric truncation radius 
is likely to be quite complicated compared to our simple power- 
law assumption of the surface density, and our comparison with 
the Herczeg et al. model is only a consistency check. 

OH lines. Sixteen rotational transitions of OH, mainly in 
the ground vibrational state, were reported by Najita et al. 
(2010) from their Spitzer observations. The data indicate high 
excitation states (/upper ~ 15-30), with equivalent energies 
E u /k up to 23,000 K above ground. OH mid-infrared thermal 
emission from our model is mainly from the hot, dense (n ~ 
10 12 -10 13 cm -3 , Tgas ~ 700-1000 K) gas very close to the 
central star at r < 0.1 AU. These lines are, however, weak 
in the 10-19 pm region and below the sensitivity of the Spitzer 
IRS instrument. OH thermal emission is also produced from 
gas in the mid-disk and outer disk regions, but this gas is 
typically at ~300 K and incapable of reproducing the observed 
high excitation lines. As suggested by Najita et al. (2010) for 
TW Hya, UV photodissociation of water is a possible source of 
the OH emission (see also Tappe et al. 2008 who suggest this 
mechanism for HH214). 

From our disk chemical network solution, we can estimate 
the total luminosity in the OH lines produced by the photodisso- 
ciation of water in the disk. We calculate the total OH line lumi- 
nosity as follows. For simplicity, in this discussion we assume 
that Lya photons ( E = 10.2 eV) are mainly responsible for the 
photodissociation of water molecules — a reasonable assump- 
tion because Lya contributes to a significant fraction of the total 
FUV luminosity (e.g., Herczeg et al. 2002; Bergin et al. 2003). 
Harich et al. (2000) studied the photodissociation of water by 
Lya photons and the photodissociation energy is determined to 
be 5.1 eV. The excess energy for a Lya photodissociation is 
therefore 5.1 eV. About 66% of the photodissociations of water 
result in OH and H, while the remaining 34% result in O and 
H atoms. Harich et al. (2000) further find that most of the OH 
products are extremely rotationally excited with a peak at a ro- 
tational level / = 45 ( E u /k ~ 45,000 K) and that almost 75% 
of the available energy goes into pure rotational excitation. This 
translates to an energy of ~4 eV per photodissociation (of H2O 
to OH) that is then eventually radiated away in a rotational cas- 
cade from the high J states (also see Mordaunt et al. 1994; van 
Harrevelt & van Hemert 2000). We calculate the energy in the 
rotational cascade in our model disk by counting the total num- 
ber of water photodissociations by photons with energy equal to 
or greater than Lya that lead to OH and multiplying this by the 
energy available for excitation (4 eV). If /VpdHLO) is the rate 
of the total number of photodissociations of water in the disk by 
Lya (and more energetic FUV) photons, then the luminosity in 
the OH cascade Loh.pd is the product 0.66 x A'pih^O) x 4 eV. 
Equating this to the observed OH luminosity then requires that 
there are ~ 1 0 4 ° photodissociations per second or that only ~ 1 % 
of the stellar Lya luminosity of 10 42 photons s -1 be intercepted 
by water in the disk (also see Bethell & Bergin 2009). If the 
solid angle subtended by the disk is ~0.1, the desired Lya flux 
absorption by water is feasible if water abundances are high 
enough (at the 10% level) to compete with the other dominant 
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sources of opacity such as dust, photoionization ionization of 
Mg, Fe, and Si, and photodissociation of CF and OH. There- 
fore, photodissociation of water can plausibly account for the 
observed line luminosities. 

In the above approximate analytic calculation, we simply used 
the Ly« flux. However, in the model we calculate the OH lumi- 
nosity due to water photodissociation using photodissociation 
rates obtained at each spatial grid cell. We use the full FUV 
spectrum, including Lya , but only count those photodissocia- 
tions where the FUV photons are energetic enough to lead to 
OH in highly excited rotational states as described above. The 
OH luminosity in the SML disk is 4 x 1 0 6 L 0 , lower than 
the total luminosity in the observed lines (upper J levels rang- 
ing from ~30 to ~ 15), 6 x 1 0 6 L 0 . We note that the cascade 
must also produce lines that are not observed in the Spitzer 
IRS band. The energy of a rotational state oc J (J + 1) and the 
cascade originates at J ~ 45. The unobserved transitions from 
J ~ 45 to J ~ 30 and J ~ 1 5 to J = 0 therefore have a 
total energy approximately 1 .5 times that in the transitions that 
are observed. The unseen luminosity must account for another 
~9 x 1 0 6 Lq, bringing the “observed” value to 1.5 x 10 -5 L 0 . 
Our SML disk results for OH produced by H 2 O photodissoci- 
ation are discrepant by a factor of ~4. Of this emission, 30% 
arises from photodissociation of water in the inner disk and the 
rest comes from the extended middle and outer disk regions to 
radii r < 40 AU (Figure 7). In the inner disk, the OH is from 
the edge of the disk near the star where the number densities are 
~10 12 -10 13 cm -3 and gas temperatures are ~500-1000 K. The 
main formation route to water is H 2 + OH and it is destroyed 
mainly by FUV photodissociation. The water abundance in these 
regions is, however low, with typical vertical column densities 
<10 15 -10 16 cm -2 in the photodissociation layer. In the mid- 
disk and outer disk regions, OH emission is dominated by the 
rim at 4 AU and beyond the rim is primarily produced from re- 
gions where the column density to the star ~(l-5) x 10 22 cm' 2 , 
densities are 10 7 -10 10 cm' 3 , and 7 gas ~ 200-600 K. In the 
mid-disk regions OH+H 2 is the dominant formation route to 
water, whereas in the cooler outer disk regions (r > 15 AU), 
H 3 0 + recombination leads to the formation of water molecules. 
Water abundances range from AdLO) ~ 1 0 6 to 10' 9 in these 
photodissociation layers, and typical water columns to the star 
are ~10 16 cm -2 . The latter is consistent with our earlier analy- 
sis, which suggested that if the water opacity were ~10% of the 
total opacity at FUV wavelengths, we would obtain OH lumi- 
nosities close to that observed. Acolumn density of ~10 17 cm' 2 
of water provides optical depth of unity at the FUV wavelengths 
needed to photodissociate water. 

If the OH is mainly prompt emission following photodisso- 
ciation of water, the SML disk appears to produce less water by 
a factor of ~4 than required by the observations. We consider a 
few mechanisms by which the water production rate may be en- 
hanced to match observations. Raising the gas temperature and 
hence increasing the water formation rate by the endothermic 
H 2 +OH reaction is one possibility and this could be achieved by 
increasing the unknown PAH abundance and FUV grain pho- 
toelectric heating in the disk. However, as H 2 pure rotational 
emission originates at the same spatial location where there is 
warm molecular gas, any increase in temperature overproduces 
the H 3 S(l) line flux. We discuss the H 2 S(l) and S(2) line fluxes 
below. 

Another possibility is that the SMH disk is the better solution 
for the inner disk mass. This disk is a reasonable match to the 
CO lines (Figure 8) and does not violate other observational 


constraints such as the lack of mid-IR water emission. Higher 
densities in the inner disk result in higher water abundances and 
the prompt OH emission increases by a factor of ~4 in the inner 
disk, bringing the total OH luminosity to ~7.8 x 1 O' 6 L 0 , and 
only a factor of two below the observed value (see Table 3). This 
model slightly overestimates the Oi 6300 A and H 2 S(2) line 
luminosities, but may nevertheless be a feasible model for the 
TWHya disk. The SMH disk also has the advantage that it better 
fits the [0 1 ] 6300 A linewidth and provides a more typical value 
for the turbulence parameter a in the inner disk (see Table 3). 

Water formation on dust grains is another mechanism that 
could increase water in the disk. Our disk models do not include 
freeze-out and desorption processes to properly consider grain- 
surface reactions. We estimate an upper limit to the effect 
of water formation on dust in a simple manner as follows. 
We assume that oxygen atoms stick on grains are instantly 
hydrogenated and released into the gas phase as water, at a 
rate equal to the collision rate of O with grains. We applied 
dust formation of gas phase water only in regions where the 
dust is cold enough to prevent thermal evaporation of O or OH 
before H 2 O is formed on the grain surface, 7 gr > 100 K. We 
use this relatively high grain temperature for thermal desorption 
so that we obtain an upper limit on the production of water on 
grains. We also require that Ay to the star is less than ~3 so that 
photodesorption can clear the water ice from the grain surface 
(Hollenbach et al. 2009). Glassgold et al. (2009) considered a 
similar mechanism to increase water production in disks but 
found that it was not significant in increasing column densities 
of warm, observable water. However, the gas temperature is 
not relevant for our calculation as the OH lines are a result of 
photodissociation and not thermal emission. Water formation of 
grains in the SML disk (with I O' 5 A7 0 within 4 AU) leads to 
a higher OH luminosity, 7 x 1 0 fl L 0 , but still a factor of ~2 
below what is observed. We emphasize that this 7 x 10 () L 0 
is an upper limit, given the optimistic assumptions on water 
formation on grains. 

Radial transport of solids or water-bearing ices in the disk, 
a process not considered by us, has been suggested as a 
possible agent causing local enhancements and depletion of 
water abundances in the disk (Ciesla & Cuzzi 2006). Water ice 
from the midplane may also be transported to greater heights 
where it is subject to thermal desorption or photodesorption, 
leading to higher water abundances in the disk surface where it 
can be photodissociated. These effects, while likely to operate 
in disks, require time-dependent models that consider radial 
and vertical transport and are beyond the scope of this paper. 
However, one simple calculation is illuminating. If water is 
transported as water ice radially inward to inner regions where 
it is photodissociated once, the mass influx needed is 10 40 H 2 O 
molecules per second. Assuming an H/H 2 O ratio of 2000 (if 
every O atom is in H 2 O ice) this translates to 2 x 10 43 H atoms 
per second or 5 x 10 7 M Q yr 1 . Since TW Hya has an age of 
10 Myr, and with a constant influx, this requires an improbable 
initial mass reservoir of 5 M 0 in the outer disk and also implies 
that the current disk will only last 10 5 years. This is unlikely. 
On the other hand, if turbulent mixing were to act vertically, 
then this has the advantage that the water can be reformed at 
the midplane and brought to the surface to be photodissociated 
repeatedly. For our adopted gas mass of 0.06 M 0 , we therefore 
require that this mass be brought to the surface of the disk 
where it is photodissociated every 10 5 years. At a typical radius, 
r = 10 AU, the distance from midplane to the Ay = 1 surface 
layer where water is photodissociated is ~ 1 AU. Thus, the 
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average vertical transport speed is only 1 AU/10 5 yr = 5 cm s _1 . 
This suggests that vertical mixing may supply the needed water. 

We have discovered a new mechanism for the rotational 
excitation of OH that may be the most promising, but difficult 
to quantify precisely at this time. Here, we do not need a greater 
production of gas phase water. The process is initiated by the 
photodissociation of OH, which leads to atomic oxygen in an 
excited electronic state (O'D) approximately 50% of the time 
(van Dishoeck & Dalgarno 1984b). The l D state of oxygen, 
due to its empty 2 p orbital, is more electrophilic than the 
triplet ground state 3 P and is highly reactive to readily undergo 
bond-forming addition reactions. O('Z)) + Hi re-forms OH very 
efficiently (with a rate coefficient yo = 3.0 x 10~ 10 cm 3 s _1 ; 
NIST Chemical Kinetics database) leading to an effective 
reduction in the overall destruction of OH. This route has no 
thermal barrier as in the reaction of the ground state 0( 3 P) 
atoms with Hi and is especially important at high gas densities 
when the rate of formation of OH by this route is faster than 
the radiative decay of the 0 1 D atom to the ground state (with 
A ~ 8 x 10~ 3 s -1 ). The fraction /q of 0 1 D that reacts with 
H2 before radiatively decaying to the ground state is given 
by fo — 1/(1 + (n cr it/«(H2))), where the critical density 
M cnt ~ yo /<! = 3 x 10 7 cm -3 . An important consequence of 
the O('D) + Hi reaction may be that this reaction produces OH 
in a rotationally excited state (with an efficiency €oh ~ 0.2 
at J ~ 20-30 states; Lin & Guo 2008) which could then 
cascade to lower J to produce the observed MIR emission by 
Spitzer. In regions with n(H2) > n cril , photodissociation of 
OH itself (rather than H2O as considered earlier) may result 
in the re-formation of rotationally excited OH. We estimate 
the contribution of this mechanism to the OH line emission as 
follows. 

Non photodissociations per second of OH will produce 
~0.5 eoh./o A'oh OH molecules in the J ~ 20-30 state per 
second. For our SML disk, we estimate that this mechanism 
could result in ~5 x 10 6 L 0 (nearly all from the inner 
r < 1 AU region), which recovers most of the observed OH 
emission which the model was deficient in when we considered 
the photodissociation of water alone. We will explore this 
mechanism for producing OH emission in greater detail in a 
future study where we will treat 0 1 D as a separate species 
for better accuracy and treat the efficiencies of eon(v, J ) 
individually. 

We conclude with a footnote to the above discussion. Mandell 
et al. (2008) recently reported the discovery of OH from two 
HAeBe stars and concluded that the emission was a result of OH 
fluorescence. The rotational energy diagram of the transitions 
observed (in this case from the v = 1 state) indicates a single 
rotational temperature which leads the authors to conclude that 
it is unlikely to be the result of “prompt” emission, i.e., OH 
produced by the photodissociation of water. Their models of 
fluorescent emission from OH match the observed data from 
HAeBe stars quite well. In their models, OH is excited to 
the v = 1 state through rovibrational transitions pumped by 
near-infrared photons and electronic transitions pumped by UV 
photons. However, the Spitzer data for TW Hya indicate mainly 
a downward rotational cascade from the high J states of i> = 0 
state (Najita et al. 2010). Therefore, in TW Hya the observed 
OH is more likely to be prompt emission. 

Hi rotational lines. Pure rotational 5( 2) and 5(1) line emission 
from H 2 was reported by Najita et al. (2010) at 4.8 x 10~ 7 
and 10~ 6 Lq, respectively. Most of the H2 emission in our 
models comes from within radii of 20-30 AU, with ~35% 


of the H2 5(2) emission and 10% of the H2 5(1) emission 
from the inner SML disk (see Figure 7). In the mid-disk 
region, Hi emission originates in the superheated dust later at 
Ay ~ 0.3-0. 1 and where the gas temperatures are 100-500 K. 
H2 formation in the mid and outer disk (r > 4 AU) is mainly 
on grains and it is destroyed by FUV photodissociation. We 
obtain H2 5(2) and 5(1) line luminosities to be 3.4 x 10 7 and 
1 .2 x I 0“ 6 L 0 , respectively. These are in reasonable agreement 
with the SpitzerlRS measurements. Bitner et al. (2008), in 
their ground-based survey for H2 using TEXES, did not detect 
H2 5(1 ) emission and obtained an upper limit of 5.0 x 10 " 7 l q . 
This suggests that half of the 5(1) emission observed by Spitzer 
IRS must originate from a more extended region outside the 
narrower slitwidth of TEXES, i.e., at r > 20 AU. Although the 
H2 5(1) emitting region is quite extended (Figure 7), half the 
H2 5(1) emission in the SML disk in fact originates at r < 
14 AU. The 5(1) flux within 20 AU is 8 x l(L 7 L 0 and the 
radial extent of emission is in moderate agreement with the 
TEXES and Spitzer data. The line luminosity in the H2 5(0) 
transition is predicted from the SML disk to be 3 x 10 " 7 Lq. 

[Oi] 6300 A and 5577 A emission. We show below that 
the [Oi] 6300 A line from TW Hya originates from the disk, 
mainly from the photodissociation of OH (also see Storzer & 
Hollenbach 1998; Acke et al. 2005) with a small contribution 
from thermal emission. The line luminosities for the [Oi] 
6300 A and [Oi] 5577 A lines from TW Hya are 10~ 5 L Q and 
1 .4 x 1 0 6 Lq, respectively (Alencar & Batalha 2002; Pascucci 
et al. 2011, S. Edwards 2010, private communication). We first 
show that the observed [Oi] 6300 A/[Oi] 5577 A line ratio of 
~7 makes it unlikely that the origin is thermal. 

The oxygen atom has a triplet ground state ( 3 P) and singlet 
D and S states which are at E u /k ~ 22,850 K and 48,660 K, 
respectively. The 0 1 6300 A transition from the 1 1) to the ground 
state and the Oi 5577 A transition from the *5 state to the 
1 D state therefore require high gas temperatures for collisional 
excitation to these high energy levels and subsequent radiative 
de-excitation. If the gas is collision dominated (in LTE), then 
the ratio of the two lines can be easily obtained from the ratio 
of n u AAE, where n u is the upper level population, A is the 
Einstein A-value, and A E is the energy of the photon emitted. 
In LTE, 77 ( 1 5) ~ (1/5 )n( l D)e AE ^ kT , where (1/5) is the ratio 
of the statistical weights and A E/k ~ 26,000 K is the energy 
difference between the two levels. Therefore 

l 6300 _ n { 1 D)A 6 iQO /S.E 6 ' iQO m -3 26000/r m 

T5577 «( 1 5)A5577AL5577 

where Ag3oo ~ 6 x 10~ 3 s -1 and A5577 ~ 1.3 s” 1 . To obtain 
the observed line ratio of 7, gas temperatures therefore need 
to be ~3600 K, which is possible with X-ray heating of the 
disk surface. However, the gas densities have to be high enough 
to produce LTE conditions. Collision partners can either be 
electrons or neutral hydrogen atoms and typical collision rate 
coefficients are y t ~5x 10 _9 (L/ 10 4 ) 06 cm 3 s -1 for electrons 
for both the transitions. The collisional rate for 0( 1 D) with 
H atoms is yn ~ 10 12 cm 3 s 1 (Krems et al. 2006). yn for 
Of 1 5) is not known. However, the rates for the relaxation of 
the O f 1 .S') state are expected to be smaller than those for the 
relaxation of the l D state, because the spin-orbit couplings in 
electronically excited atoms are generally weaker (R. Krems 
2011, private communication). We conservatively assume yn ~ 
10- 12 cm 3 s -1 for Of 1 5) as well. Therefore, electrons dominate 
collisions if x e > 4 x 10 4 . For levels to achieve LTE, electron 
densities n e > A/y e ~ 10 8 cm~ 3 are needed for the 5577 A line 
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Figure 9. Contours of the thermal line ratio 0 1 6300 A/0 1 5577 A as a function 
of the gas temperature and electron density from a statistical equilibrium 
calculation of the O atomic levels are shown. TW Hya has an observed ratio 
of 7 (solid line) which requires high electron densities at typical neutral gas 
temperatures. For an electron abundance x e ~ 0.01 and T ~ 8000 K, gas 
densities n ~ 10 9 cm -3 are needed to reproduce the observed 6300 A/5577 A 
line ratio by thermal O i emission. 

( n e > 10 6 for the 6300 A line). Collisions with H atoms require 
neutral densities about 2500 times these values. Even for gas 
with high electron fractions x e ~ 0.1, gas densities therefore 
need to be high, «h // 1 0 9 cm' 3 . These conditions are likely to 
be attained only in the inner disk, where gas at the surface may 
be dense, be partly ionized with sufficient neutral oxygen, and 
at T > 3600 K. 

For gas at lower densities a balance between radiative de- 
excitation and collisional excitation and de-excitation of the 
different levels determines the level populations. We calculate 
the line ratios as a function of electron density and gas temper- 
ature and Figure 9 shows the contours for the observed value 
of 7 and other values about this number. From the contour plot, 
one can see that even for relatively high gas temperatures of 
~8000 K the required electron densities are ~10 7 cm -3 im- 
plying gas densities of ~10 9 cm -3 for an electron fraction of 
x e ~ 0.01, which is typical. Such conditions are unlikely, ex- 
cept perhaps in the innermost regions, r < 0. 1 AU. We also note 
that soft X-rays that can heat the gas to such high temperatures 
typically have very small absorption columns ~ 1 0 1 9 cm' 3 , and 
very little mass is expected in gas at these high temperatures 
and densities. Such low masses may not be able to reproduce 
the observed high optical O i line luminosities. 

For typical densities in the hot disk surface, n ~ 
10 5 -10 7 cm -3 , and electron fractions x e ~ 0.0 1 — 1 0 3 the elec- 
tron densities are n e < 10 5 cm' 3 , and the observed line ratio 
requires T > 30,000 K. We consider this extremely unlikely as 
it is very difficult to heat disk gas to these high temperatures 
and moreover retain a significant amount of neutral oxygen. 
The recent calculations of Ercolano & Owen (2010) also show 
that the 5577 A line is appreciable relative to the 6300 A line 
only when gas densities are high, as in their primordial disk 
models where the Oi emission is from the inner r < 1 AU 
region. In their models, the O i emission is thermal in origin and 


they quote Z.6300/T5577 ~ 13-50 for the primordial disk models 
and ~300-700 for the transition disk models, supported by the 
above analysis. Only one of their many models, a primordial 
disk with L x ~ 2 x 10 3 ° erg s' 1 has a low line ratio ~ 13, where 
the emission presumably originates from gas with high electron 
density. 

In the SML disk, thermal emission contributes to about 20% 
of the observed [Oi] 6300 A line and 10% of the 5577 A 
line. In the inner disk, the luminosities are 1.2 x 10' 6 Lq 
and 10 " 7 a 0 for the thermal component of the 6300 A and 
5577 A lines, respectively, and the emission comes from regions 
where n e ~ 10 5 — 10 8 cm' 3 and T ~ 4000-8500 K at the disk 
surface. In the SMH disk, the calculated thermal [0 1] 6300 A 
line luminosity is higher, 2 x 10' 6 L 0 (Table 3). There is also 
thermal [Oi] 6300 A emission from the mid and outer disk 
surface that is heated by X-rays to T > 5000 K, the luminosity 
here is ~ 1 . 1 x 1 O' 6 L 0 . There is negligible thermal [0 1] 5577 A 
emission from the mid and outer disk regions. 

We propose here that most of the 0 1 forbidden line emission 
from TW Hya is the result of photodissociation of OH in 
the disk. The OH photodissociation layer in the SML disk 
is almost spatially coexistent with the H2O photodissociation 
layer, and this generally is located near the H/H2 transition in 
the disk. OH photodissociates about 50% of the time to produce 
an oxygen atom in the 1 D state and occasionally (5% of the 
time) in the higher '5 state (e.g., Festou & Feldman 1981; 
van Dishoeck & Dalgarno 1984a). Subsequent electronic de- 
excitation to the ground state produces the oxygen doublet at 
6300 A and 6363 A ( l D— 3 P) and the 1 .S’— 1 D transition produces 
the 5577 A emission. The atoms in the O l D state produce the 
6300 A+6363 A doublet in the branching ratio of 3 (Storey & 
Zeippen 2000). The ratio of these line strengths is the same as the 
branching ratio because both lines originate from the same upper 
state. These calculations are supported by comet observations of 
the photodissociation products due to solar Ly a (Morgenthaler 
et al. 2007). We again count the number of photodissociations 
of OH in the model disk as we did for water, and multiply this 
by the energy of the transition and the theoretically calculated 
branching ratio. 

Our inclusion of the 0(*D) + H2 — »■ OH reaction leads to a 
lower effective photodissociation rate of OH as discussed earlier. 
If this reaction were not considered, then prompt emission 
produces a 6300 A/5577 A line ratio ~7. However, in regions 
where the densities are high, such as in the inner disk, the 
line ratio is expected to be lower as the reactive nature of 
the 0 1 D atom will result in a suppression of 6300 A emission. 
Although the 1 S state of oxygen is also more electrophilic than 
the triplet ground state, the unpaired electrons likely make the 
Of 1 ,S')+H2 reaction less favorable. Further, the radiative decay 
of O'S is faster (A5577 > 200 A 63 oo) requiring even higher 
densities for the reaction with H2 to be competitive. The 5577 A 
line is therefore not expected to be lowered significantly. In 
dense regions of disks, the ratio of 6300 A/5577 A is therefore 
expected to be lower than ~7. /V ()H photodissociations per 
second of OH will produce ~0.5(1 — /oWoh 6300 A photons, 
where (1 — /o) is the reduction factor due to 0( 1 D) + H2 —*■ 
OH which suppresses the radiative transition leading to the 
6300 A line. We therefore calculate that the [Oi] 6300 A and 
[Oi] 5577 A line luminosities due to OH photodissociation are 
1 . 1 x 10~ 6 L q and 9.3 x 10” 7 Lq in the inner disk, respectively, 
and these lines originate in the same region as the OH cascade, at 
/• < 0.1 AU, where n ~ 10 12 -10 13 cm' 3 and T ~ 500-1000 K. 
OH photodissociation in the mid-disk regions produces line 
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luminosities 5 x lO" 6 L 0 (6300 A) and 7 x 10" 7 L 0 (5577 A). 
This emission is produced from gas with n ~ 10 7 — 10 10 cm -3 
and T ~ 200-600 K out to radii r ~ 40 AU. Therefore, most of 
the O i luminosity from TW Hya comes from photodissociation 
of OH. Ercolano & Owen (2010) recently propose that the Oi 
forbidden line emission seen in disks may arise from an X-ray 
driven wind. Our picture differs from theirs in two ways — in our 
model the O i emission arises through a non-thermal process and 
the emission arises from the disk atmosphere rather than in a 
wind. This perspective is consistent with the Oi 6300/5577 A 
line ratio and the lack of a blueshift for the emission. 

The total luminosities from the SML disk (thermal +photodis- 
sociation) are therefore 8.4 x 10 -6 L f;) for the [Oi] 6300 A line 
and 1.6 x 10“ 6 L 0 for the [Oi] 5577 A line. Our 6300/5577 A 
line ratio is ~6. The SMH disk which was a better match for the 
OH emission gives slightly higher [O i] emission, 1.5x10 " 5 L q 
for the 6300 A line and 1.92 x 10 6 Lq for the 5577 A line for 
a ratio of 7.8 (also see Table 3). 

The intrinsic FWHM linewidth of the observed O i lines are 
of order 10 km s -1 (Pascucci et al. 2011). This corresponds 
to a velocity dispersion of about 5 km s -1 . In the SML disk, 
the Oi lines arise mostly in the mid-disk regions where the 
projected Keplerian speeds are low, <2 km s . In addition, 
the gas is relatively cool and turbulent speeds are expected to 
be low (Hughes et al. 2011). Hence, the contribution to velocity 
dispersion from thermal and turbulent velocity fields is expected 
to also be low, <1 km s 1 . However, the photodissociation 
of OH leads to a recoil of the resultant 0( 1 D) or 0( 1 5) 
atoms at speeds of order 2 km s -1 , and the [Oi] photons are 
emitted prior to collisional relaxation (van Dishoeck & Dalgarno 
1984a). We estimate the linewidth of the Oi 6300 A line by 
adding the various components and the Keplerian, thermal and 
recoil contributions as appropriate. We obtain a linewidth of 
~6 km s -1 , narrower that the observed value. In the SMH disk, 
the calculated linewidth is more in agreement with observations, 
~11 km s -1 , because of the increased prompt emission close 
to the star, although the total Oi 6300 A line luminosity is 
overproduced by a factor of ~ 1.5 (Table 3). 

We end with a caveat that although the above production 
mechanism via OH photodissociation explains the O i 6300 A 
and O i 5577 A line strengths, ratio, and lack of an observed 
blueshift in TW Hya, it may not be a more general explanation 
valid for all disks in which optical forbidden line emission 
has been detected, especially if these lines are blueshifted and 
originate in a flow. For disks with different surface density 
distributions and FUV and X-ray luminosities, there may be a 
more dominant thermal contribution to [O i] emission. Further, 
our disk models are static and do not include any emission 
that might arise in an extended wind (e.g., Ercolano & Owen 
2010). Although we expect that our disk density and temperature 
structure is reasonably accurate within our grid, we do not 
account for any emission that might arise from heights above 
our vertical extent of z — r. Densities in these regions are 
typically low, and usually not a significant source of emission. 
For the optical lines (O i and S ii, discussed below) however, the 
upper states are at ~20, 000—40, 000 K and the emission is very 
sensitive to temperature, and the contribution from the wind in 
this case might be significant. 

S ii optical lines. Pascucci et al. (2011) re-analyze the optical 
data of Alencar & Batalha (2002) and report the detection of 
the [S ii] 4069 A line (~3 x 10“ 6 L 0 ) from TW Hya. They 
do not detect other lines of S ii and place an upper limit of 
1.7 x 10- 7 Lq on the 6731 A line. We find that in our SML 


disk, most of the S ii emission arises from the inner disk in 
hot, dense gas that also results in the thermal Oi emission 
calculated above. Our calculated Sii line luminosities are 
3.9 x 1O” 6 L 0 for the 4069 A lines, and 8.6 x 10- 7 ,5.6x 10~ 8 , 
and 1.1 x 10- 7 L 0 for the undetected lines at 4076 A, 6716 A, 
and 6731 A, respectively. Our [Sii] 4069 A line is slightly 
stronger than observed, however, sulfur may be somewhat 
depleted compared to our assumed abundance (2.8 x 1 0' 5 ). 

Nen and Nem. [Nen] emission from TW Hya is strong 
and has been detected by Ratzka et al. (2007), Herczeg et al. 
(2007), Pascucci & Sterzik (2009), and Najita et al. (2010). 
There is some discrepancy between the measured values, with 
more flux observed by the Spitzer study. Najita et al. (2010) 
conclude that [Nen] emission arising from a more extended 
region ( r > 10 AU) explains their higher flux at all three 
observational epochs, which the narrow slits of the ground- 
based observations exclude (Pascucci & Sterzik 2009). A revised 
calibration (Pascucci et al. 2011) recovers the Spitzer flux and 
suggests that most of the [Ne ii] emission comes from within 
~10 AU. Earlier work by Herczeg et al. (2007) measured a 
line luminosity of 4.4 x 1 0 6 L 0 with their ground-based study, 
similar to the Spitzer data. There is therefore some indication 
that the [Nen] line flux (or the underlying continuum) might 
be variable (Pascucci et al. 2011). Here we use the Spitzer IRS 
measurements of Najita et al. (2010) for the [Nen] and [Nem] 
emission, (3.9 — 4.7) x 10 -6 L 0 and 2 x 10~ 7 Lq. 

[Nen] emission from the disk models is found to be quite 
insensitive to the surface density distribution. This is because in 
X-ray heated regions, the [Nen] emission mainly depends only 
on the intercepted X-ray photon flux and the temperature of the 
emitting gas (HG09). A simple estimate of the expected [Nen] 
line luminosity can be made using the analytical expression 
derived in HG09, 
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where (px ~ 10 39 s 1 is the 1 keV X-ray photon luminosity. 
In the above, we also assume that the disk intercepts a fraction 
/ ~ 0.7 of the stellar X-ray flux (our grid in the numerical 
work typically extends to z max = r). A typical gas temperature 
in the [Ne ii] emitting layer is T ~ 2000 K. The expected [Ne ii] 
luminosity for the TW Hya disk is therefore 4 x 1 O' 6 Lq. 

From our full, detailed numerical calculations, we obtain 
total [Ne ii] and [Ne hi] line luminosities for the SML disk as 
3 x 10~ 6 Lq and 3 x 10 - 7 Lq , respectively. Approximately 25% 
of the [Nen] comes from gas in the inner disk (see Figure 7). 
The mid-disk region (3.5 AU < r < 20 AU) contributes the 
remaining 75% to the [Nen] luminosity in the SML disk. The 
[Ne ii] arises from the X-ray heated and partially ionized surface 
of the disk where gas temperatures range from 1000 to 4000 K. 
[Nen] emission from the SML disk is restricted radially to 
r < 10 AU (Figure 7) and peaks strongly at the vertically 
extended 4 AU rim. The X-ray-produced [Nen] is somewhat 
lower than the observed line flux, which could be explained 
either by the presence of an EUV component and/or by 
variability of the stellar high-energy photon flux. We note that 
the accretion rate of TW Hya is observed to be variable on 
timescales similar to that for the [Nen] line flux (Eisner et al. 
2010), ~a few years. As the X-ray and any [Nen] producing 
EUV radiation fields may be a result of accretion-induced stellar 
activity (Robrade & Schmitt 2006; HG09), this may naturally 
explain the variability of the [Nen] line. 
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Pascucci & Sterzik (2009) resolve [Ne ii] in their spectra, and 
the observed linewidths and blueshifted emission are consistent 
with a photoevaporative flow from the star. The line profile is 
very well reproduced by EUV ( hv > 13.6 eV) photon-driven 
photoevaporation models (Alexander 2008). However, we do 
not include EUV in our source spectrum for TW Hya and can 
recover most of the observed flux from X-ray heated gas alone. 
The remaining ~25% of the observed Nen flux may be due to 
the EUV layer. We can thus estimate the EUV photon luminosity 
as being <3 x 10 40 s _1 . The [Nen]/[Nein] line ratio may 
provide some measure of the origin of the emitting gas (the 
completely ionized EUV-heated layer or the partially ionized 
X-ray-heated layer below the EUV-heated layer; 10,000 K in the 
EUV case and ~ 1000-3000 K for X-rays), but is not definitive 
because of the uncertain nature of the EUV spectrum (HG09). 
However, the most natural explanation is an origin in the X-ray 
heated gas, where rapid charge exchange of Ne ++ with H atoms 
quenches the [Nem] line and produces the observed [Nen]/ 
[Nem] ratios. (For a more detailed discussion, see HG09.) 
The observed blueshift may also arise for an X-ray induced 
photoevaporative flow (Gorti et al. 2009; Owen et al. 2010), 
and the lack of need for significant EUV to reproduce the 
[Ne n] suggests that any photoevaporative flows may originate 
in primarily neutral, X-ray heated gas. 

The far-infrared and submillimeter emission lines discussed 
hereon originate only in the mid and outer-disk regions and are 
the same for the two disk models, SML and SMH. The line 
luminosities are listed in Table 3. 

[Oi ] and [ C n] fine structure lines. We obtain a [O i] 63 /im 
line emission of 6 x 1 0 6 L 0 , a factor of ~2 higher than the 
Herschel PACS measurement of Thi et al. (2010). Our calculated 
[Oi] 145 /im line luminosity is 3 x 10 ~ 7 L q , and consistent 
with the upper limit of 4 x 10 " 7 ^o obtained by Thi et al. 
(2010). Emission is extended and produced by ~ 5 0-1 00 K gas 
at r ~ 30-120 AU. 

Our model overproduces [O i] 63 jx m emission, and we 
suggest that inclusion of water freezing on grains might alleviate 
this problem. Freezing of water ice on grains colder than 100 K 
removes gas phase elemental oxygen from the outer disk where 
dust grains are cold. Photodesorption would release some of this 
ice and at Ay < 3, all water ice is photodesorbed from grains 
(Hollenbach et al. 2009). Including these processes reduces the 
[O i] 63 /tm emission from the disk compared to a disk with 
no ice formation, to 4 x 10~ 6 Lq, in better agreement with 
the PACS observations. The [O i] 145 jrm line is reduced from 
3 x 10~ 6 Lq to 2 x 1 0 6 Lq. A more detailed treatment of 
freeze-out of molecules on dust grains and desorption processes 
will be a subject of future work. Another possibility is that 
the PAH abundance in the outer disk is even lower than we 
assume, with a depletion factor of nearly ~ 1 000. This lowers 
the temperature in the outer disk sufficiently and decreases the 
[Oi] 63 jim line luminosity to 4 x 1 0 6 Lq. If freezing were 
also included, this would bring the O i emission even closer to 
the observed value. However, we note that decreasing the PAH 
abundance everywhere in the disk alters the vertical structure 
of the mid-disk region (r < 20 AU) sufficiently to deteriorate 
fits to observation of other species such as H 2 , OH, and Oi 
forbidden lines. There may, in fact, be a radial gradient in the 
PAH abundance in the disk, with a higher degree of depletion in 
the outer disk, if PAHs were to enter the disk frozen out on dust 
grains as suggested by Geers et al. (2009). Another possibility is 
that the PAH abundance is low everywhere in the disk, but that 
other sources of heating such as mechanical viscous dissipation 


are effective at the surface in the disk(r < 30 AU, e.g., Glassgold 
etal. 2004). 

Thi et al. (2010) modeled [Oi] emission from TW Hya 
and concluded that the observed flux together with their 13 CO 
observations indicate a low gas disk mass by factors of 10-100 
compared with that expected from the disk dust mass and 
assuming a gas /dust ratio of 1 00. However, we find that lowering 
the gas disk mass by 10 and retaining the dust mass only 
marginally improves model fits to the observed [Oi] emission, 
and worsens the other line emission fits considerably. Most 
of these discrepancies can be attributed to the increased grain 
absorption cross section per H in a disk with a low gas/dust 
ratio. H 2 emission increases (FUV heating by PAHs becomes 
significant) by a factor of ~2-3, OH prompt emission is reduced 
due to the increased (relative to gas) absorption of Lyor photons 
by dust, and there is excess emission from CO pure rotational 
lines. The HCO + 4-3 line requires high densities for emission, 
and this line luminosity is almost a factor of 10 lower for the 
low mass disk than the observed value. Increasing the PAH 
abundance by 10 as in the preferred model of Thi et al. would 
raise gas temperatures even higher and further increase H 2 and 
CO emission. Although we do not concur with the Thi et al. 
(2010) results, we note that our calculated model line fluxes are 
in good agreement with results from other models they consider. 
Our model is most similar to their low PAH abundance models 
(their series 2 and 3) and if they were to base their conclusions 
on these alternate models, a disk mass close to ~0.06 M 0 is 
inferred from their results. For these series, they also seem to 
find that the 0 1 line luminosity is quite insensitive to disk mass 
in the range we consider. Their conclusion of a low gas disk mass 
appears to be largely based on the ratio of the 12 CO/ 13 CO line 
fluxes. In order to accurately interpret this measurement, a more 
accurate disk chemical model which includes 13 CO as a species 
is necessary. Our model does not include 13 CO. We also note that 
Thi et al. assume a constant 13 CO to 12 CO ratio to calculate the 
13 CO flux and do not include a proper chemical model of 13 CO 
as well. Moreover, the 13 CO/ 12 CO line ratio will be significantly 
affected by processes such as selective photodissociation (Lyons 
& Young 2005) and freeze-out processes, which are not included 
in our disk models or those of Thi et al. (2010). It remains to 
be seen whether a low gas disk mass is needed to reproduce the 
observed CO isotope flux ratio. 

The calculated [Cii] 158 jim line luminosity is 3 x 10 - 7 Lq 
and is consistent with the upper limit of 4.5 x 10 7 Lq set by 
the Herschel observations. 

CO rotational lines. The observed line luminosities of the 
CO 6-5, 3-2, and 2-1 lines are 5.0 x 10” 8 ,2.5 x 10 -8 , and 
7.9 x 10~ 9 Lq, respectively, whereas we obtain 8 x 10~ 8 , 2.4 x 
10 -8 , and 7.9 X 1O- 9 L 0 for these lines. Purely thermal line 
broadening in the outer disk provides better correspondence 
between data and model results. This is in fair agreement with 
the results of Qi et al. (2006) who assume only a small amount 
of turbulence (AU ~ 0.3 km s~ 3 ) in their analysis and with 
the more recent results of Hughes et al. (2011) who find that 
the turbulent linewidths are low. A power-law surface density 
profile also matches the data, but has to be arbitrarily truncated 
at ~70 AU (GH08). We consider it more realistic to assume that 
the surface density declines more smoothly with radius, as in a 
disk evolving viscously and subject to photoevaporation (also 
see Hughes et al. 2008; Isella et al. 2009). 

CO rotational line emission constrains the distribution of 
gas in the outer disk. We do not model disk masses higher 
than O.O6M 0 , because we expect that appreciably higher disk 
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masses will make the disk gravitationally unstable. Lower disk 
masses by a factor of 10 ( M<ji S k = 0.006 M 0 ) do not fit the 
CO rotational line emission well. We find that although it is 
possible to match the observed CO 2-1 and 3-2 line emission 
by altering the surface density profile in the outer disk for the 
lower mass disk, the CO 6-5 line emission is higher by a factor 
of ~4 than the SML value, and hence discrepant with data. We 
expect freezing will not affect the emission line luminosities 
calculated. A simple estimate made by discounting contributions 
from spatial locations where the dust temperature is <20 K 
only slightly lowers the calculated line luminosities by ~15%. 
Photodesorption may also keep the CO molecules from freezing 
on dust grains in the surface regions of the outer disk. We note 
that we do not need severe radial depletion of CO to explain the 
observed line ratios (e.g., Qi et al. 2006), although some degree 
of depletion (less than a factor of three) in the mid-disk region 
may help lower the CO 6-5 line luminosity from the model disk. 

HCO + 4-3 line, van Zadelhoff et al. (2001) measured the 
line intensity of the HCO + 4-3 transition from TW Hya using 
the James Clerk Maxwell Telescope (JCMT), and the line 
luminosity is 7.7 x 10“ 9 L 0 . Our calculation yields 4 x 10“ 9 L 0 , 
nearly a factor of two lower than observed. We could obtain a 
better match by changing the surface density distribution in the 
outer disk to increase the strength of the HCO + 4-3 line, but this 
then also changes the CO rotational line intensities and makes 
that agreement poorer. Further, increasing the outer disk surface 
density also results in a high overall disk mass, bringing the 
disk closer to being gravitationally unstable. The HCO + line 
mainly originates in relatively warm (T ~ 20-100 K) and high 
density gas (n > 10 7-8 cm -3 ) at r < 50 AU. The abundances 
of HCO + calculated by the model in the emitting layers of the 
disk are a factor of 20-70 lower than the typical dark cloud 
abundance of 5 x 10“ 9 adopted by van Zadelhoff et al. (2001) in 
their analysis, and thus consistent with their derived “depletion” 
factors of 10-100. 

We conclude this section with the caveat that our models may 
not be unique solutions to the gas surface density distribution in 
the disk. There are theoretical uncertainties that are inherent 
in chemical rates and the microphysics, as well as several 
important processes (such as freezing of ices) have been either 
simplified or ignored. The complicated and computationally 
intensive modeling procedure also makes a full exploration of 
parameter space impossible at this time. Nevertheless, we find 
it promising that theoretical models with reasonable parameters 
can reproduce observed emission ranging from the optical to 
the submillimeter wavelengths with an accuracy to within a 
factor <2. 

4. THE ORIGIN OF THE INNER GAS HOLE 


dust grains into larger dust particles and planetesimals. Grain 
growth is expected to be more efficient in the inner disk as the 
densities there are higher (e.g., Lissauer & Stevenson 2007). 
The population of small dust decreases due to grain growth, and 
as larger solids are not efficient at absorbing starlight because of 
their lower collective surface area, this process is manifested as a 
decrease in the dust continuum emission at shorter wavelengths 
where the warm small dust in the inner disk would emit. Grain 
growth does not remove disk gas, and no depletion in gas surface 
density within the inner regions is expected (unless gas giant 
planets have already formed, which is the next scenario). In this 
case, there would be continued accretion onto the star (as is 
observed for TW Hya). 

Our modeling of the gas disk around TW Hya when compared 
with observational data does not support the grain growth 
explanation alone as a mechanism to cause the dust hole. Based 
on the constraints set by gas line and continuum emission, we 
argue that gas must also be depleted in the inner disk region 
(r < 4 AU), with a surface density decrease ~l-2 orders of 
magnitude compared to what one would obtain by extrapolating 
the surface density from the outer optically thick, massive disk. 
Although grain growth may explain the decrease in dust opacity 
inside 4 AU, a mechanism that depletes gas surface density by 
~10-100 is also necessary. 

Giant planet formation. A widely favored explanation for 
the observed dust depletion in transition disks is the creation 
of a gap by a giant planet, followed by the formation of a 
hole as the inner disk accretes onto the star. A gas giant planet 
exerts dynamical torques on the disk and clears material in its 
immediate vicinity to create a gap (e.g., Artymowicz & Lubow 
1996; Bate et al. 2003; Rice et al. 2003). The disk interior to the 
planet is expected to rapidly accrete onto the central star. The 
outer disk mainly accretes onto the planet, with a small fraction 
of the mass flow (depending on the ratio of planet mass to disk 
mass) continuing to accrete onto the star (Lubow & D’Angelo 
2006). Pressure maxima in the disk just outside the planet may 
act to filter dust particles (Rice et al. 2003) and further deplete 
the inner disk of dust. Our modeled surface density distribution 
supports this explanation for the TW Hya disk, with the potential 
planet located between 2 and 4 AU where the surface density 
may be quite low (with gas in accretion streams) and consistent 
with the presence of a gap. The accretion streams then feed a 
depleted inner accretion disk at <2 AU. 

We can use the surface density distribution determined by the 
gas emission models to set constraints on the mass of the planet 
orbiting TW Hya. As the planet must open a gap in the disk, or 
at least establish a large surface density contrast at its location, 
we can use the gap-opening criterion to obtain a minimum mass 
M p of the planet (e.g., Lin & Papaloizou 1986) 


Modeling of gas emission lines and comparisons with ob- 
servational data indicate that the surface density of gas is also 
depleted in the inner r < 4 AU dust depleted region. The models 
are further consistent with a massive, optically thick outer disk 
beyond the 4 AU rim. At least three different explanations have 
been proposed to explain the dust holes in transition disks, of 
which TW Hya is a leading example (e.g., Najita et al. 2007). 
The distribution of gas in the disk that we derive may help 
distinguish between these various scenarios. We examine the 
implications of each of the explanations for transition disks and 
discuss our model results in each context. 

Grain growth with no gas depletion in the dust hole. The first 
step toward planet formation in disks via the core accretion 
scenario is the coagulation of primordial sub-micron-sized 
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where H is the scale height of the outer disk at the radius 
r ~ 4 AU. From our disk models, H/r ~ 0.04 at 4 AU, and 
therefore the planet has a minimum mass of ~0.03-0.3 Mj for 
a ~ 10“ 3 to 10“' in the outer disk. Simultaneously, one can 
appeal to recent hydrodynamical models of accretion onto and 
past a massive planet (Lubow & D’Angelo 2006), which find 
that a planet with at least this minimum mass must accrete matter 
at a rate ~10 times higher than the rate onto the star. Given the 
current age of the system (~ 10 Myr), the planet is likely to have 
undergone accretion for at least 1 Myr. Even with the lower 
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estimate of the stellar accretion rate ~4 x 10~ 10 M 0 yr -1 , this 
implies a minimum planet mass M p ~ 4 Mj . An upper limit to 
the planet mass can also be estimated from the stellar accretion 
rate. The fact that gas in the outer disk can still flow through the 
gap induced by the planet suggests that the planet is not massive 
enough to form a tidal barrier and repel disk matter. Depending 
on a, this limits the planet mass to ~7 Mj{a = 4 x 10~ 3 ) or 
~2 Mj(a — 4 x 10 4 ) (see D’Angelo et al. 2010). We conclude 
that a plausible planet mass is likely to be ~4-7 Mj. 

The mass of the planet inferred above is high enough that 
dust particles with sizes larger than ~1 /iin may be trapped in 
the pressure maximum in the disk at the outer edge of the gap 
(Rice et al. 2003). Dust particles with sizes of ~2 /i m have been 
modeled to exist in the inner disk (C02), consistent with this 
scenario. Since the inferred dust distribution in the outer disk 
has a maximum size of a max = 1 cm, and only sub-micron-sized 
grains cross the gap, dust filtering would result in a decrease in 
dust mass inside the gap. Therefore an increase in the gas/dust 
mass ratio across the gap by a factor of ~ 100 (since the dust 
mass is proportional to V( a max)) is expected, which is consistent 
with our SMH model disk. 

A potential issue with the planet hypothesis is the migration 
rate. For this case, we need to consider Type II migration as the 
planet must be massive enough to create a gap. Viscous diffusion 
of the disk causes the planet to drift toward the star efficiently, 
but only as long as the local disk mass ( nr 2 'L(r ), r ~ 4 AU) is 
greater than the planet mass. For the mid and outer disk surface 
density distribution, the local disk mass is ~1 Mj. Thus, the 
expected planet of ~4 My is stable against migration. 

Photoevaporation. The strong FUV and X-ray flux from 
TW Hya should drive massive photoevaporative flows from the 
disk (Ercolano et al. 2008; Gorti & Hollenbach 2009; Gorti et al. 
2009; Owen et al. 2010). Photoevaporating flows have been ob- 
served by resolved ground-based [Ne ii] 12.8 // m emission from 
TW Hya (Pascucci & Sterzik 2009). The observed blueshifted 
line profiles are well reproduced by theory (Alexander 2008; 
Ercolano & Owen 2010). 

However, whether the inner hole in TW Hya is carved by 
photoevaporation remains a question. The location of the radius 
of the hole ~4 AU is consistent with theory as it is larger than 
the so-called critical radius, ~ 1 AU where the hole first appears 
in photoevaporative theory before the disk is eroded radially 
outward by direct rim illumination. However, photoevaporation 
theory traditionally predicts that a “clean” hole, completely 
devoid of dust and gas, must rapidly form on viscous timescales 
(e.g., Alexander et al. 2006). The viscous clearing timescale 
in the inner disk is f v ; s ~ r 2 /v or t vis ~ */(GM*r)/(ac 2 ). For 
the TW Hya disk, f v i S ~ 10 4 years at 4 AU for a ~ 0.01. 
Therefore, the surface density drop of ~ 100 at 4 AU implies 
an e-folding factor of five, or ~5 x 10 4 years after gap opening 
by photoevaporation. For a higher surface density in the inner 
disk and a drop of ~10 at 4 AU, gap opening must be even 
more recent. Although there must be ongoing photoevaporation 
from the disk as is observed, our results appear to refute current 
photoevaporation theories for the creation of the inner hole in 
TW Hya, unless the gap creation epoch is extremely recent, 
which seems unlikely given the ~ 1 0 7 year age of TW Hya. 

We note that one possibility where one could avoid a 
special observing epoch and where photoevaporation could 
work as a gap creation mechanism would be if the thermal 
photoevaporative flow that emanates from the rim and disk 
surface is partially re-captured by stellar gravity to form the 
inner disk and eventually accrete onto the star. In our models. 


the photoevaporative mass loss rate is ~4 x 10” 9 M 0 yr -1 . 
For a stellar accretion rate of ~4 x 1O~ 1O M 0 yr -1 , 10% of 
the flow needs to be recaptured by the star. We use simple 
timescale arguments to show that this is a feasible process. For 
recapture, the flow timescale it/) to the sonic radius r s where gas 
escapes needs to be comparable to the viscous timescale f v ; s of 
the photoevaporating gas. As tf ~ r s /c s , and using r s — r g / 2 
and the expression for ? v j s in the preceding paragraph, we obtain 
tf/ty i S ~ air g /r) l l 2 . For FUV/X-ray driven photoevaporation 
where the flow is launched subsonically (GH09, GDH09) 
from the rim at r = 4 AU, and a typical gas temperature 
of ~2000 K, r g = 31 AU. For 10% of the flow is to be 
recaptured due to viscosity, we need tf/t v - IS ~ 0.1 or a value 
of the viscous parameter a = 0.1. If the accretion rate is 
~1O” 9 M 0 yr -1 , the mass captured needs to be ~25% for 
the calculated photoevaporation rate. Therefore, for this re- 
capture process to be possible, a needs to be relatively high 
in the inner disk, ~0. 1-0.25. However, this process is a 
complicated phenomenon that requires full two-dimensional 
hydrodynamical calculations to resolve, including turbulent 
viscosity in the flow that can redistribute the angular momentum. 
This problem will be the subject of future study. 

We also note that even if the inner opacity hole is not created 
by photoevaporation, the process must still act to disperse the 
disk. The observed [Ne n] line profiles indicate this is an ongoing 
process in the TW Hya disk. We can estimate the lifetime 
of the remaining disk around TW Hya using FUV/X-ray 
photoevaporation rates from theory (GH09) and the results from 
the SML disk. From our calculated disk photoevaporation rate 
4 x 1CD 9 M q yr -1 , we derive a disk lifetime of M^/M pe as 
~10 7 years. Note that with the reduced disk mass of Thi et al., 
the expected disk lifetime is < 1 0 fi years, which places us in 
an unlikely viewing epoch and therefore seems implausible. We 
calculated a time evolution model of the disk, with the SML disk 
as an initial configuration and assuming a viscosity parameter 
a = 0.005 in the outer disk using the photoevaporation models 
described in Gorti et al. (2009). When the effects of viscous 
evolution are included, the remaining lifetime of the disk is 
calculated to be 5 x 10 6 years. We note that a higher assumed 
value of a would lead to even shorter disk lifetimes. 

Most of the photoevaporation occurs at the 4 AU rim, 
although the mass loss rate derived is lower compared to 
that from holed disks with directly illuminated rims (e.g., 
Gorti et al. 2009; Owen et al. 2010) because of the partial 
shielding of the rim in this case by gas in the inner disk. 
The photoevaporation rate is higher than recent measurements 
of the stellar accretion onto TW Hya, which range from 
4 x 10~ 10 to 3 x 10~ 9 Mq yr -1 (Muzerolle et al. 2000; Alencar 
& Batalha 2002; Eisner et al. 2010). This is consistent with 
the requirement for photoevaporation to create a gap in the 
disk, when the accretion rate at ~1— 4 AU must drop below the 
photoevaporative mass loss rate. 

If photoevaporation has caused the inner hole and not a giant 
planet, then the amount of gas we infer to be present in the inner 
disk and the observed accretion onto the star imply that either 
the gap has been very recently created (< 1 0 5 years ago) or that 
the gap was created earlier and that the observed accretion is 
caused by viscosity in the flow causing significant re-capture of 
the photoevaporating gas and accretion onto the star. If a giant 
planet has caused the inner hole, then the accretion rate in the 
outer disk is expected to be at least a factor of 10 higher than 
M*, and in this case the derived photoevaporative mass loss rate 
is lower than the outer accretion rate, which is consistent with 
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the giant planet hypothesis and, in this case, not consistent with 
the pure photoevaporation scenario for hole formation. 

5. CONSTRAINTS ON THE DISK MASS FROM LINE 
EMISSION MODELING 

In this paper, we model the gas disk of TW Hya and conclude 
that the best model which reproduces the observed gas emission 
line luminosities is one where there is an optically thick disk 
(0.06 Mq) extending from —3.5 AU outward, with a depleted 
inner region that contains — 1 0 4 to 10“ 5 M 0 of gas extending 
to the star. The modeling suggests a lower gas mass in the inner 
disk by ~ 10-100 than obtained by extrapolating the surface 
density of the outer disk radially inward. Most of the disk mass 
(0.06 Mq) is contained within 100 AU, and the surface density 
in the outer disk declines sharply beyond this radius. The gas/ 
dust mass ratio is consistent with the canonical ratio of 100 in the 
outer optically thick disk and is enhanced by factors of —5-50 
in the inner disk where both gas and dust are severely depleted 
compared with the extrapolated surface density of the outer 
disk. Although the outer disk gas/dust ratio is the same as the 
interstellar ratio, the cross-sectional area of dust per H nucleus 
is — 100 times less than interstellar dust, because the disk dust 
size distribution extends to much larger sizes, a max ~ 1 mm, 
than interstellar dust where a max — 0.2 /xm. Therefore there are 
far fewer small dust particles ( which provide most of the cross- 
sectional area) per H nucleus, since much of the dust mass is 
tied up in the largest size (~1 mm) particles. 

Inner disk. The observed CO rovibrational emission and limits 
on HtO line emission between 10 and 19 /xm constrain the inner 
disk mass and are best explained by a gas mass — 10“ 5 M G 
within r < 3.5 AU. Lower disk masses (1O“ 6 M 0 ) are very 
discrepant with observed high J CO rovibrational emission and 
can be ruled out. A full, undepleted disk with mass — 10“ 3 M G 
is rejected based on the expected continuum emission due to gas 
opacity and the lack of observed mid-IR water emission. A disk 
mass of 10“ 4 M 0 is, however, also in reasonable agreement with 
data. There is some uncertainty in the inner disk dust mass as 
inferred from observations (~2 x 10“ 8 Mq, C02; — 10“ 10 Mq, 
Eisner et al. 2006), but the conclusion that the hole is depleted 
of gas and dust compared to the outer disk appears quite robust. 

The amount of mass present in the inner disk is fairly 
consistent with the observed accretion rate. Assuming steady 
state accretion rate onto the star AL = 3 it vX, where the 
viscosity v = etc 2 /LI. a being the viscosity parameter, c s 
the sound speed of gas, and Q. is the angular velocity (e.g., 
Pringle 1981), we can estimate a from the observed M* and 
from X and c s derived from our models. For an accretion rate 
M* = 4 x 10“ 10 M q yr” 1 , and at r = 1 AU, we then obtain 
a ~ 0.04 for the SML disk and -0.004 for the SMH disk. 
(If M* were —10 “ 9 M 0 yr -1 , then a — 0.1 for the SML 
disk and a — 0.01 for the SMH disk.) At the inner edge, 
close to the magnetospheric truncation radius of 0.06 AU, 
a is higher, —0.08-0.008. Our derived values for a for the 
preferred SML disk are somewhat higher than the 0.01-0.005 
value typically used for protostellar disks, but consistent with 
the higher gas temperatures and ionization levels expected in the 
inner disk(King et al. 2007; Turner et al. 2010). For the SMH 
inner disk mass of 10“ 4 Mq, the inferred value of a is more in 
line with values typically used. 

Outer disk. HCO + and CO rotational lines place upper limits 
on the amount of gas that is present in the outer disk and 
are consistent with a disk gas mass —0.06 Mq. All the other 
observed line emission originates mainly in the disk at radii 


r < 30 AU, which contains about 30% of the total disk mass. 
We find that decreasing the gas mass by a factor of 10 lowers the 
HCO + luminosity by a factor of five due to its origin in dense 
gas. Since our model already underestimates this luminosity 
by a factor of two, this would produce a major discrepancy in 
the model luminosity versus the observed HCO + luminosity. In 
the adopted surface distribution with a tapering outer edge, CO 
pure rotational lines are also quite sensitive to the mass in the 
outer disk. The lower mass disk results in a high CO 6-5/2-1 
line ratio of 16, compared with the observed ratio of 5.8 (see 
Table 1). Increasing the gas mass by a factor of 10 brings the 
disk mass to nearly the stellar mass, and such a disk would 
be highly gravitationally stable (Pringle 1981). We therefore 
consider a more massive disk than 0.06 Mq very unlikely. In the 
outer disk, gas line emission may therefore constrain disk gas 
mass to within a factor of < 10, even though most of the gas line 
emission originates in the surface layers which do not trace the 
bulk of the disk mass. 

We can also set limits on the disk mass from arguments based 
on the presence of a planet and from expected mass loss rates due 
to photoevaporation. If a planet were present and accreting at 
— 10 times the stellar accretion rate, this would imply a relatively 
short disk survival time for a 0.006 Mq disk. Even for a low 
accretion rate onto the star of M ; , = 4 x 1O“ 1O M 0 yr -1 , the 
disk survival time would be —6 x 10“ 3 M 0 /(1OM») — 1 Myr, 
which is short compared to the age of the system. This situation 
is made even more dire if the disk mass were lower or if the stellar 
accretion rate is at the higher range of measurements. Further, if 
the disk is losing mass at the estimated photoevaporation rate of 
4 x 10“ 9 Mq yr -1 , a 0.006 M 0 disk would dissipate in <1 Myr, 
which seems unlikely given the age of TW Hya. In view of all 
the above, our preferred gas disk mass is 0.06 Mq. 

6. CONCLUSIONS 

We infer the gas surface density distribution of the transition 
disk around TW Hya by modeling gas line emission from 
the disk, including H 2 , HCO + , and CO rotational lines, CO 
rovibrational lines, OH line emission, and [Neil], [Nem], [Oi] 
fine structure emission and optical 0 1 and S 11 forbidden line 
emission. Our standard best-fit model (SML) disk reproduces 
observed line emission well to within a factor of —2 for most 
lines, except OH emission which is discrepant by a factor of 
—4. We consider prompt emission following photodissociation 
of water as mechanism to produce OH. If we can increase our 
gas phase water production at the disk surface by a factor of a 
few, we can bring all emission lines in the model to within a 
factor of two of the observations. Vertical advection of icy grains 
in the mid disk region to surface layers where the ice mantles are 
photodesorbed may be the most likely source of the extra water 
vapor. We further speculate that photodissociation of OH itself 
may result in re-formation to a rotationally excited state in very 
dense regions, which might lead to the observed OH emission. 
The alternate model of a higher inner disk mass (SMH) provides 
good matches to the observed emission as well. The SMH disk 
is better at reproducing the OH emission (prompt) and the [0 1 ] 
6300 A linewidth, although the CO rovibrational line emission 
is a poorer fit and some emission lines are overproduced. Our 
main results on gas line emission are as follows. 

1. The model disk that best matches the observational data 
indicates a significant decrease in gas surface density 
in the inner opacity hole that causes a dip in the mid- 
infrared continuum. Our best models have a gas mass 
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of 1CT 4 to 10” 5 M o at 0.06 AU < r < 3.5 AU and 
0.06 M 0 between 3.5 AU < r < 200 AU. We exclude a 
full undepleted disk with 10~ 3 M Q and inner disk masses 
of 10 6 M e and below. The gas/dust mass ratio in the inner 
disk is inferred to be greater by a factor of ~5-50 than the 
canonical value of 100 which applies in the outer disk. 

2. CO rovibrational emission and upper limits on fUO rota- 
tional lines provide good constraints on the mass in the inner 
disk. The observed CO rovibrational emission is modeled 
to arise from hot (T ~ 500-1000 K) gas inside ~0.5 AU, 
from disk surface regions where hydrogen is predominantly 
atomic. 

3. The mass of the outer disk is constrained to lie within a 
factor of a few of 0.06 M 0 . Significantly higher masses 
lead to gravitational instability. Significantly lower masses 
lead to very poor fits to the CO rotational spectrum and to 
the HCO + J = 4-3 line luminosity. In addition, given the 
relatively rapid accretion and photoevaporation of the outer 
disk, lower disk masses give unreasonably short outer disk 
lifetimes (<1 Myr). 

4. We propose that OH mid-infrared emission originates 
from prompt emission following photodissociation of water 
in the disk very close to the star ( r < 0.1 AU) and 
from the mid-disk regions ( r ~ 4-40 AU). Our model 
underestimates the production of OH by a factor of four, 
and we speculate that vertical mixing may bring more water 
to the surface to be photodissociated. Grain formation of 
water ice followed by photodesorption could also increase 
the water vapor at the surface and thereby also increase the 
OH emission. However, enhanced water production might 
also lead to enhanced [O i] 6300 A forbidden line emission 
by a factor of ~3 above the observed value. We suggest 
another promising mechanism for producing OH emission 
in dense regions where the photodissociation of ground 
state OH leads to an excited 0 1 1) atom that reacts with H 2 
to form rotationally excited OH that cascades giving rise to 
the observed MIR lines. 

5. We argue that the forbidden lines of [Oi] at 6300 A and 
5577 A arise due to the photodissociation of OH mainly 
from the mid-disk regions ( r ~ 4-20 AU). Such an origin 
can successfully explain the Oi 6300 A/5577 A ratio of 
~7 observed. About ~20% of the [Oi] 6300 A emission 
is thermal in nature and has equal contributions from the 
inner disk (r < 4 AU) and the outer disk (r ~ 4-10 AU). 

6. Ne 11 emission may not require the presence of a significant 
EUV flux from TW Hya and can be explained by X-ray 
heating of gas. If EUV is present, we estimate a photon 
luminosity of <3 x 10 40 s -1 in order to explain the slight 
deficit in X-ray-produced Ne 11 line fluxes. 

7. H 2 pure rotational line emission is extended and is mainly 
from radii r ~ 4-30 AU. [0 1 ] 63 /zm emission, CO pure 
rotational lines, and HCO + 4-3 line trace the outer disk at 
r ~ 30-120 AU. 

We conclude that grain growth alone cannot be responsible 
for the dust hole in the TW Hya disk because gas line emis- 
sion modeling indicates a factor of ~ 10-100 depletion in the 
gas surface density as well. The transition disk morphology of 
TW Hya may be explained by the presence of a ~4-7 M j planet 
at r ~ 2-3 AU that causes the change in surface density distri- 
bution near this region. The disk is presently photoevaporating 
with an estimated mass loss rate of 4 x 10 -9 M 0 yr _I . If there 
were no planet, then photoevaporation alone may be in the pro- 
cess of developing a hole. The calculated photoevaporation rate 


is greater than the stellar accretion rate, suggesting that photoe- 
vaporation could have created the gap. However ongoing stellar 
accretion at current rates poses problems for present photoevap- 
oration theories. Regardless of the hole-creation mechanism, 
the remaining disk lifetime due to photoevaporation and vis- 
cous evolution is estimated to be ~5 million years. 
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